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§ 1. InTRODUCTION 

THe term epitaxy (‘arrangement on’) was introduced by Royer (1925) 
to denote the phenomenon of the oriented growth of one crystal upon 
another. The first examples of such growth were observed to have 
occurred naturally on minerals; two crystals of different species grow 
together with some definite and unique orientation relationship between 
their crystal axes. Many such cases of oriented intergrowths have been 
observed, and comprehensive lists of the examples have been given by 
Wallerant (1902) and by Mugge (1903). 

Such observations led to experiments designed to bring about oriented 
growth in the laboratory. The first successful attempt appears to be 
that of Frankenheim (1836), the previous claims of Wakkernagel (1825) 
having been disputed. Frankenheim found that sodium nitrate grew 
from solution in parallel orientation on the surface of a calcite crystal. 
and also found that potassium iodide, bromide and chloride were oriented 
when crystallized from solution on to a mica cleavage surface. The first 
systematic experiments in which a series of related structures are grown 
upon each other were made by Barker (1906, 1907, 1908), who studied, 
amongst other things, the growth of alkali halides upon each other. 
Barker found that some alkali halides oriented upon each other, while 
others did not. It was noticed that oriented growth is more likely to 
occur if the molecular volumes of the two alkali halides are nearly equal. 
Royer (1928) repeated much of the work of Barker, and extended it to 
cover many other examples of growth from solution. In the meantime. 
the discovery of x-ray diffraction had greatly increased knowledge of 
the structure of crystals. On the basis of his results, Royer put forward 
three rules of epitaxy, the most important of which is that oriented 
growth occurs only when it involves the parallelism of two lattice planes 
which have networks of identical or quasi-identical form and of closely 
similar spacings. The difference between the network spacings is usually 
expressed in terms of the percentage misfit, which is defined as 100(b—a)/a. 
a and b are the corresponding network spacings in the substrate and 
overgrowth respectively. The experiments indicated that the misfits 
should be no more than about 15%, this limit being most strikingly 
illustrated by the growth of alkali halides both upon alkali halides and 
upon a mica cleavage surface. 

The discovery of electron diffraction at about the time when Royer 
carried out his experiments greatly increased the scope of studies of 
epitaxy. Whereas previous work had. to be carried out by means of 
optical microscopy, with its inherent limitation on the types of layers 
which can be studied, the electron diffraction technique provided a means 
of studying very thin surface films. It was now possible to examine 
deposits other than those grown from solution, and to study earlier stages 
of growth. Although the early development was slow, it soon became 
apparent that epitaxy occurs with a wide variety of deposits on single 
crystal surfaces. Very thin chemically grown layers (Aminoff and Broome 


Study of Epitaxy in Thin Surface Films 175 


1936, 1938, Wilman 1940), electrodeposited metal films (Cochrane 1936, 
Finch and Sun 1936) and metal layers condensed from the vapour phase 
(Bruck 1936) are all found to be oriented in some cases. One of the 
most striking features of the results of these new experiments is that 
Royer’s rule concerning a good numerical fit between the substrate and 
the overgrowth is not always obeyed. This applied particularly to the 
case of metals condensed on to heated rocksalt surfaces (Bruck 1936), 
for which misfits as big as 38% were observed. Attempts were made to 
explain the occurrence of these large misfits by consideration of the 
formation of initial layers which are oriented with a low misfit, and upon 
which other orientations formed at a later stage. A mechanism of this 
kind, in which the occurrence of twinning played a vital role, was 
considered in detail by Menzer (1938 a, b, c) for the case of silver and 
nickel on rocksalt. 

A further concept which was introduced by Finch and Quarrell (1933, 
1934) as a result of electron diffraction studies is that of ° basal plane 
pseudomorphism ’. This involves the formation of an initial oriented 
film which has an abnormal crystal structure. The bulk structure of 
the deposit material is constrained so that the lattice plane parallel to 
the surface is identical in size to that of substrate. Frank and van der 
Merwe (1949 a, b, c) followed up this concept, and developed a theory 
of epitaxy which involves, as a necessary condition of orientation, the 
formation of an initial pseudomorphic layer. This theory predicts 
limiting misfits of the same order of magnitude as those observed by 
Royer. However, certain difficulties arise concerning the stability of 
such layers (Smollet and Blackman 1951), and more recent experiments 
have failed to confirm the necessity of a small misfit (e.g. Schulz 1951 b, 
1952 a). In addition, these more recent experiments have included the 
examination of extremely thin layers (mean thickness less than 104), 
and this has provided much needed evidence of the structure of the 
initial nuclei of an oriented deposit. No support for the concept of 
basal plane pseudomorphism is found. 

The main purpose of this article is to give an account of the electron 
diffraction evidence concerning the formation of an oriented overgrowth, 
and to consider the relation of this evidence to existing theories. This 
will include a discussion of the uses and limitations of electron diffraction 
as a technique in such studies. It has not been possible to include a 
detailed review of the results obtained by other methods, but such 
results have been included and discussed where they are directly 
related to the electron diffraction work. Comprehensive reviews 
containing most of the results obtained by optical microscopy, together 
with the few results obtained by other methods such as x-ray diffraction 
will be found elsewhere (Neuhaus 1950-51, 1952, van der Merwe 
1949, Seifert 1953). Any consideration of the growth of deposits 
on amorphous substrates, upon which fibrous orientations are often 
observed has also been omitted. 
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§ 2. THe Meruops or StupyING Epitaxy 
2.1. The Optical Microscope Method 


The optical microscope has been used almost exclusively for the study 
of growth from solution. The essential requirements of the method are 
that the crystallites of the deposit should be sufficiently big, and should 
have a well-defined external form so that the nature of the boundary 
faces can be determined. The orientation of each crystallite is then 
determined separately, and any preference can be established. 


\ 


2.2. The X-ray Diffraction Method 


The x-ray method involves a direct determination of the orientation 
of the overgrowth crystallites, independent of their external form. 
A composite diffraction pattern from a large number of crystallites is 
obtained, and any preferred orientation can be directly established. 
Various techniques are possible, but it is usually necessary to use a 
rotation method, so that any orientation of the overgrowth will give 
rise to observable patterns. Almost any type of layer can be examined 
in this way, provided that it has sufficient thickness. 


2.3. The Electron Diffraction Method 


A great deal of work has been carried out during the last twenty-five 
years on the application of electron diffraction to the study of epitaxy. 
Most of the comments in this paper refer to the use of fast electrons 
(30 to 100 kev), most of the work having been carried out in this energy 
range. 

The most common and more generally useful method is the reflection 
method, in which the overgrowth is examined in situ on the substrate. 
The transmission method may also be used; there are two distinct 
techniques. Firstly, a single crystal substrate in the form of a thin 
film may be used, and the overgrowth formed on this film (see, for example, 
Wilman (1940) or Finch and Sun (1936)). Secondly, the overgrowth is 
formed on a single crystal substrate and then detached (e.g. by dissolving 
away the substrate) so as to form a thin film which can be examined by 
transmission (see, for example, Bruck (1936) or Shirai (1937)). 

In most cases, the electron beam illuminates a specimen area of at 
least 0-lmmX0-lmm, so that the diffraction pattern arises from a 
large number of overgrowth crystallites. In some electron diffraction 
instruments, particularly those combined with electron microscopes. much 
finer and more intense electron beams are used. The number of crystallites 
irradiated may then be much smaller, and in certain cases of high resolution 
examination it would be possible to study the individual crystallites. This 
technique has not yet been applied to study of epitaxy. 

Electron diffraction has one considerable advantage over other methods. 
Even quite thin surface films (see § 2.3.6) give very intense patterns, 
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which are clearly visible on a fluorescent screen. A thorough visual 
examination of the overgrowth can therefore be made at various crystal 
settings, and photographs can be readily obtained at all appropriate 
settings. 

Further details of the general technique of electron diffraction may be 
found elsewhere (Beeching 1936, Finch and Wilman 1937, Pinsker 1953, 
Raether 1951, Thomson and Cochrane 1939). 


2.3.1. Types of Specimen Suitable for Examination 


In the case of transmission specimens, the thickness of the films must 
not be more than a few hundred angstréms. This limit arises because 
of the inelastic scattering processes, which give rise to a continuous 
background on the diffraction patterns. Since studies of epitaxy necessitate 
the use of single crystal substrates, this thickness condition provides a very 
severe limitation. There are two possible ways of increasing the maximum 
specimen thickness. Firstly, the energy of the electrons can be increased. 
According to Mollenstedt (1946), a tenfold increase in film thickness is 
possible, for polycrystalline aluminium, if the energy is increased from 
80 to 500 kev. Secondly, the inelastically scattered electrons may be 
filtered from the diffraction pattern. Boersch (1953) has shown that this 
results in a considerable improvement in the contrast of the pattern, but 
no data are available concerning the application of the technique to the 
examination of thick films. 

For reflection work, the substrate should have a macroscopically flat 
surface, preferably parallel to a definite lattice plane. Cleavage faces of 
a crystal are very convenient, particularly since they usually provide 
surfaces which are fairly flat on an atomic scale, at least over a limited 
area. It is often necessary, particularly in the case of metal substrates, 
to use surfaces prepared by some form of etching. Only in very rare 
cases does the etching process produce surfaces as flat as a cleavage face ; 
generally the surface is rough on an atomic scale. The suitability of 
such a surface as a substrate depends upon various factors such as the 
extent to which the contact planes may be determined (see §§ 2.3.2 and 
2.3/3). 

The overgrowths may be prepared by any method which provides 
erystals which give sufficiently bright diffraction patterns. 


2.3.2. Examination of the Surface Form of the Substrate 


Any theoretical discussion of the mode of orientation of an overgrowth 
must involve a consideration of how the two lattices are joined together. 
Hence it is important to establish the planes of contact between the 
two lattices. It has been shown by several experimenters (see, for 
example, Pashley (1952 a) that the relative orientation does depend upon 
these contact planes. The first step in their determination involves the 
investigation of the topography of the substrate. 
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The ideal type of substrate to use is one which has an atomically 
smooth surface. The contact planes can then be deduced immediately 
from the relative orientation. In practice, the nearest approach to such 
a surface is a cleavage face of a good single crystal. These usually have 
a step structure, the regions between the steps consisting of surface 
which is fairly flat on an atomic scale. Observations for mica are given 
by Tolansky (1948, p. 122) and for sodium chloride by Schulz (1952 a). 

A substrate which cannot be prepared in the form of a cleavage surface 
is usually quite rough on an atomic scale. It is important that the 
topography of such a surface be examined. The only methods of much 
value are electron microscopy or electron diffraction, although optical 
techniques such as the Bridgman (1925) method can be used to identify 
well-developed crystal facets. Electron microscopy has not so far been 
used very much in this field. This is partly because the limit of resolution 
of replica techniques has been, until recently, no better than about 1004. 
There now exists the carbon replica technique which has been developed 
by Bradley (1954), and which is claimed to have a much higher resolution. 
It is to be hoped that more extensive use of electron microscopy will 
follow, although contrast limits might well prove to be serious. 
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Refraction displacement as a function of the grazing angle of incidence (a). 
=1-0002 ; camera length 40 em. 


Electron diffraction is capable of supplying a certain amount of 
information about surface topography, which can be quite detailed in 
favourable cases. The degree of flatness of a surface can be judged 
from the amount of displacement caused by the refraction effect. For 
fast electrons (30 to 100 kv) the refractive indices of materials varies 
approximately from 1:00005 to 1-:0004. Because of this low value the 
refractive displacement is measurable only if the electron beam enters 
the surface at a small glancing angle x. The variation of the displacement 
with « is given in fig. 1, for a value of »=1-0002. The displacement is 
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expressed both in degrees and in terms of the shift obtained on a photo- 
graphic plate at a typical distance of 40 cm from the specimen. In normal 
circumstances this shift can be detected only if it is greater (at best) than 
about 0-2mm. The value of » is usually no more than about 5°, so that 
a measurable refractive displacement is observed if the surface is perfectly 
flat (see fig. 2(a)). If the surface is rough on an atomic scale, the 
electron beam enters the actual surface at a large glancing angle, and 
no appreciable refractive displacement results (see fig. 2 (c)). An inter- 
mediate effect occurs with a surface which is fairly smooth, but which 
contains shallow undulations (see fig. 2 (b)). If the mean slope of the 
undulations is less than about 5°, then observable displacements occur. 
For a parallel incident electron beam, the glancing angle of entry « will 
vary according to the part of the surface entered ; a range of refractive 
displacements occurs, causing the diffraction spots to be tailed down 


Refraction effects and surface form. 
(a) Atomically flat surface ; (b) wavy surface ; (c) rough surface. 


towards the shadow edge of the surface. Such patterns are commonly 
found from electrolytically polished surfaces (Kranert, Leise and Raether 
1944). The refractive displacement therefore gives useful information 
concerning surface flatness. Some confusion occurs in the literature 
concerning the meaning of the terms ‘ atomically smooth ’ and ‘ atomically 
flat’. Surfaces, such as that of fig. 2 (a), which show the full refraction 
effect are usually considered to be flat to the extent that their boundary 
is made up of an almost completely filled atomic plane, although the 
extent to which this need be so is unknown. In this article the term 
‘atomically flat’ will only be used to denote surfaces which fulfil this 
condition of full refraction. Surfaces which give rise to tailing (fig. 2 (b)) 
will not be included in this category, although some authors refer to 
these and similar surfaces as ‘ atomically smooth’. Some authors also 
deduce that surfaces are ‘ atomically smooth’ from the fact that sharp 
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Kikuchi lines are observed. This conclusion is, in the writer’s opinion, 
completely unjustified because the formation of the Kikuchi lines has 
no direct relationship to the surface form. 

The electron diffraction technique cannot be used to determine the 
area of the ‘ atomically flat ’ patches (e.g. areas between cleavage steps) 
on a surface; it can only be used to provide an approximate lower 
limit to this area. The relevant condition to be fulfilled is that the 
electron beam should both enter and leave the surface through the 
atomically flat region. Since the path length of the electrons in the 
material is limited, this will certainly be fulfilled if the patches are more 
than a few hundred angstréms in extent in the direction of the electron 
beam. Thus the minimum requirement of a surface, such that it fulfils 
the definition of ‘atomically flat’ given above, is that it should be 
made up of a large number of regions which are flat enough to give rise 
to the full refraction displacement, the linear size of these regions being 
a few hundred angstréms. 

When the substrate surface is of the type shown in fig. 2 (c), it is of 
some importance to establish whether there are any well-developed 
facets, since any overgrowth will be formed upon these facets. If a 
surface is bounded by a large number of similar facets, these can be 
identified in favourable cases. There are two distinct effects. 

Firstly, if the incident beam is so directed that electrons can enter or 
leave the surface via these facets at a small glancing angle, a refractive 
displacement will occur. This will shift the diffraction spots on the 
photographic plate in a direction normal to the plane of the facet, so 
that identification is possible. Such discrete shifts are not often observed. 
Instead, diffraction spots are * tailed’ along the inwards normal to a 
particular set of facets, and such * tailing’ is usually attributed to 
refraction effects. Several difficulties arise in this interpretation (Pashley 
1951 a); in particular, some effect such as a wavy facet surface has to 
be postulated, in order to account for the observation of tailing rather 
than discrete spot displacements. 

The second effect concerns the form of the diffraction spots. The 
theory has been developed by Laue (1936, 1937 a, 1948) and Laue and 
Riewe (1936). Lf a small crystallite has well-developed boundary faces, 
the reciprocal lattice points corresponding to the crystallite will be 
extended in directions normal to those boundary faces. These extensions 
or ‘ stacheln ’ give rise to the splitting of diffraction spots, so that several 
spots appear instead of one. Each subsidiary spot arises from an inter- 
section of the Ewald sphere with the ‘ stacheln’ through the corresponding 
reciprocal lattice point. From the form and magnitude of the splitting 
the boundary facets may be identified. For a given system of facets, the 
pattern may be calculated from a formula derived by Pashley (1951 a). 

The application of this theory to observed practical cases involves 
certain interpretative difficulties. The first examples of such effects 
were for the growth of silver layers on rocksalt (Lassen 1934, Kirchner 
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and Lassen 1935) and for the growth of cobalt layers on copper (Cochrane 
1936). In both of these cases Laue (1936, 1937 a) deduced, from the 
observed diffraction patterns, that the (111) type planes were the 
prominent boundary planes. From the geometrical viewpoint this 
interpretation was satisfactory, since the observed patterns were 
consistent with the occurrence of extensions along the [111] directions 
in the reciprocal lattice. However, Kirchner and Rudiger (1937) have 
criticized the facet interpretation, and disagreement with Laue was 
expressed (Kirchner 1937, Laue 1937 b). There are three main points : 
(1) The splitting was accompanied by features which indicated the 
occurrence of repeated (111) twinning in the layers. The splitting effects 
can be associated, at least in part, with this twinning (Cochrane 1936). 
(2) Kirchner and Rudiger (1937) maintain, from other electron diffraction 
features, that the cube planes and not the (111) planes, form the boundaries. 
(3) The observed effects correspond to ‘ stacheln’ which are longer and 
more intense than those to be expected from Laue’s theory, for crystals 
of a reasonable size. 

More recently (Pashley 1951 a) splitting effects have been observed 
from etched silver crystals, for cases where twinning was definitely absent. 
The length and intensity of the ‘ stacheln ’ were more near to the theoretical 
prediction, but the ‘ stacheln ’ were found to be asymmetrical about the 
centres of the reciprocal lattice points. No satisfactory explanation for 
this has been found, although surface deformations of a special type 
might be responsible. 5 

To summarize, there exists, in theory, a method of detecting small, 
numerous surface facets by means of a study of the form of diffraction 
spots. This method usually gives rise to certain difficulties of interpre- 
tation, which make the identification of facets not completely certain. 

-More work on the interpretation of fine structure features is necessary ; 
in particular, correlation with electron microscope observations might 
prove to be helpful. 


2.3.3. The Determination of Contact Planes 


When a rough single crystal surface is used as a substrate, an overgrowth 
will often occur in several orientations, corresponding to deposition on the 
various surface facets. It is then necessary to determine which of these 
planes is giving rise to a particular overgrowth orientation. In favourable 
cases, unambiguous interpretations may be arrived at (Pashley 1952 a) 
as a result of electron diffraction examination. A careful study is made 
of the surface at different azimuth settings. When a facet is set approxi- 
mately perpendicular to the electron beam, the growth on that facet 
gives rise to a prominent secondary diffraction pattern, which results from 
scattering by both the substrate and the overgrowth. When the same 
facet. is set nearly parallel to the electron beam, only a primary pattern 
is obtained from the overgrowth. These differences depend upon 
whether or not the electron beam traverses a substrate projection in 
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entering or leaving a thin deposit on a facet surface. By correlating 
information of this kind obtained at different crystal settings, the number 
of possible contact planes can be reduced considerably. This number is 
further reduced by considering as possible contact planes only those pairs 
of lattice planes in the substrate and overgrowth which are parallel. The 
extent to which this restricts the choice depends considerably upon the 
symmetry of the two structures, and upon their relative orientation. 

By means of this method, the true contact planes have been deduced 
with certainty for silver halides on etched silver crystal surfaces (Pashley 
1952a) and for cadmium oxide on electropolished cadmium crystal 
surfaces (Lucas 1952). 


2.3.4. Spacing Measurements and their Limitations 


Since it is difficult to measure the electronic wavelength directly, 
accurate spacing measurements are often carried out with the aid of a 
calibration specimen. This technique can be successfully employed for 
both reflection and transmission specimens; provided adequate pre- 
cautions are taken, the inaccuracy of calibration need be no greater 
than about 0-1°% for reflection and can be much less for transmission. 
The choice of a standard calibration material is of importance. Thallium 
chloride prepared by deposition in vacuo is very suitable and reliable 
(Boswell 1950, 1951). 

The accuracy of measurement is, in most cases, limited by the lack of 
sharpness of the diffraction spots. This restriction becomes particularly 
serious in the reflection case, where effects due to surface form and lack 
of penetration are often considerable. Diffraction spots are broadened or 
elongated, sometimes in an asymmetrical manner, so that the point to 
which measurements should be made is ill-defined. An absolute accuracy 
of about 1°, is usually possible for spacings of planes perpendicular to the 
surface, but for the planes parallel to the surface the spacings can sometimes 
be measured only to a much poorer accuracy. Only in a comparatively few, 
although sometimes vital, cases is it impossible to distinguish between the 
spacings of the substrate and the overgrowth (further discussion is given 
in § 4.4). 


2.3.5, Continuous Examination during the Growth of a Layer 


When an oriented layer is formed in situ in the evacuated diffraction 
camera, it is possible to obtain a continuous record of the change of the 
diffraction pattern during the growth. This technique has been applied 
mainly to the study of the growth of layers from the vapour phase. Two 
distinct methods are used. In the first (Schulz 1951 b, 1952 a), a number 
of different specimens are prepared in the evacuated camera and these are 
then examined. Further deposits on the same specimens may be also 
examined, In the second (Uyeda 1942) a single specimen is used and 
deposition is carried out while the electron beam is being used for the 
examination of the surface. The advantage of this latter method is that 
all stages of growth are studied in one experiment, and that the same 
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area of one specimen is examined. Effects due to specimen patchiness 
or variation from one specimen to another, are avoided. The disadvantage 
however, is that in some cases the electron beam is known to influence the 
growth of the deposit (Collins and Heavens 1952). We have found, in 
this laboratory, that growth upon several alkali halides is seriously 
influenced by the electron beam. The effect seems to be due to some 
permanent damage to the substrate, probably arising from a partial 
decomposition at the surface. In such cases, the first method must be 
used. 

The importance of these techniques is that the growth may be carried ° 
out in a very controllable manner, so that the very early stages of an 
oriented layer may be studied. The value of the techniques depends upon 
the extent to which electron diffraction can be used in the examination 
of extremely thin layers. 


2.3.6. The Limit of Detection of a Surface Film 


Since a transmission specimen of only 100 4 in mean thickness will yield 
a bright diffraction pattern, it is to be expected that a uniform film of 
thickness ¢ on a flat substrate surface will give an observable pattern if 
the path of the beam traversing the film is at least 100A. For a grazing 
angle of incidence «, this condition becomes ¢/sin «>100A, which 
reduces, for a typical value of sin «=0-02, to t+>24. A single uniform 
atomic layer should therefore be observable. The precise limit of 
detection would depend upon the deposit material used, and would be 
a function of both the angle « and the incident electron energy H. The - 
dependence upon # is difficult to predict, since it is controlled by an 
intensity balance between the inelastic scattering from both the substrate 
and the overgrowth, which causes a general background, and the elastic 
scattering from the overgrowth. For a rough substrate, the sensitivity 
of detection is lower because the electron path length in the overgrowth 
is smaller for a uniform deposit of a given thickness (see Raether 1949, 
1950). A non-uniform film is also less easy to detect than a uniform 
layer. 

Any experimental investigation concerning these conclusions is difficult 
because it involves independent measurements of film thickness in the 
range 1 to 104. Weighing methods are too insensitive, and most other 
recognized methods are also unsuitable. Various authors have attempted 
to measure the thickness of evaporated films by methods involving an 
estimation of the amount of material which leaves the source. In several 
cases bright diffraction patterns have been reported from layers of only 
a few angstroms in thickness. Schulz (1951 b, 1952 a) has used such a 
method, to an estimated accuracy of about 20%. As a result, he was 
able to show that oriented layers of alkali halides as thin as 1 or 24 give 
observable patterns. 

This aspect has been studied in more detail (Newman and Pashley 1955) 
with the aid of a radioactive tracer technique. Atomically flat silver 
crystal surfaces were used as substrates, and two types of layer were 
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investigated : (1) radioactive metallic copper deposited from the vapour : 
(2) silver bromide formed by the chemical attack of radioactive bromine 
vapour. The f-activity of the specimens was used to give thickness 
measurements for deposits as thin as 0-14. The results showed that 
layers of mean thickness less than 1 A gave observable diffraction patterns, 
the limit of detection being about 0-84 for the copper and 0-4A for the 
silver bromide, for electrons of energy about 45 kev. Detailed qualitative 
information concerning the effect of the grazing angle of incidence « was 
obtained. 

An important feature of these latter experiments is that in neither case 
did the deposits grow as uniform layers, but in the form of separated 
crystallites distributed over the surface. This aggregation has little effect 
upon the limit of detection, provided the crystallites are sufficiently small, 
as they are in this case. If the crystallite size is such that the incident 
coherent electron intensity is considerably reduced as the beam traverses 
it, the sensitivity of detection is reduced. Since the sensitivity of detection 
is greatest for a uniform deposit, it follows from the results with copper 
and silver bromide that for substances of average atomic number, and 
for atomically flat substrate surfaces, partial as well as complete atomic 
monolayers may be detected and examined. 

In cases where growth occurs in the form of a large number of isolated 
crystallites, detection may be possible only after the crystallites have 
reached several atomic layers in thickness, although the average thickness _ 
of the deposit material may still be less than that of one atomic layer. 

~The mode of growth of such crystallites cannot be studied in detail during 
the early stages. 


2.3.7, The Form of the Overgrowth Crystallites 


Although electron diffraction can be used to provide useful information 
concerning the form (e.g. size and shape) of the overgrowth crystallites, 
great caution is necessary in the interpretation of the results because of 
the limited penetrating power of the electrons. For this reason, it is also 
difficult to make any reliable deductions, unless the overgrowth is formed 
on an atomically flat surface. Even in this latter case, the overgrowth 
must be sufficiently thin to ensure that the diffraction arises from the 
whole of the layer. This normally implies that the layer must be no more 
than about 104 in average thickness. If the thickness of the individual 
crystallites exceeds about 104A, the information will be further limited 
when the lateral extent of the crystallites exceeds about 100A. 

When layers of only a few angstréms in mean thickness are studied by 
the methods outlined in § 2.3.5, the mode of nucleation of the deposit 
‘may be investigated. If the layer spreads over the surface so as to form 
a coherent film of almost uniform thickness, the diffraction pattern 
consists of spots which are considerably elongated in a direction perpen- 
dicular to the surface. In the limiting case of a single atomic layer, this 
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elongation results in the formation of continuous streaks, corresponding 
to diffraction by a two-dimensional grating. If the layer forms as a 
number of separated crystallites, the pattern consists of distinct spots. 
The detailed structure of these spots is related to the size and shape of 
the crystallites. The interpretation of the effects is exactly similar to 
that used for the determination of surface form (see § 2.3.2). A check 
upon this kind of electron diffraction evidence has recently been made, 
in one isolated example. Miyake (1938) deduced that antimony oxide 
(Sb,03) grew on a heated stibnite surface in the form of well-defined 
prisms which were oriented in a regular manner. Ito, Ito and Watanabe 
(1954) have since studied the same tarnishing reaction by electron 
microscopy, and confirm Miyake’s result, which had been based entirely 
upon electron diffraction evidence. More work of this kind would be 
very helpful, in view of the various sources of uncertainty. 

If the crystallites of the deposit grow such that they are partially 
embedded below the surface of the substrate, the size and shape of the 
primary diffraction spots correspond to the part of the crystallite protruding 
above the surface. The embedded part might, however, contribute to the 
secondary diffraction patterns. It is possible for measured sizes of 
crystallites to be false because of crystal irregularities (e.g. stacking 
faults) which reduce the size of good coherent crystal. 


§ 3. THE OBSERVED CASES OF EPITAXY 


As a result of the many studies which have been made, a large number 
of cases of epitaxy has been observed. A comprehensive list of examples 
has been given by van der Merwe (1949) and further data will be found 
in the other reviews quoted in §1. A survey will now be given of all 
of the cases of epitaxy which have been observed by electron diffraction 
methods, and these results will be compared with the corresponding 
observations which have been carried out with other techniques, mainly 
optical microscopy. 


3.1. Alkali Halides upon Themselves 


The first extensive experiments on the growth of alkali halides upon 
each other were carried out by Barker (1906, 1907, 1908). His results 
are included in table 1, and refer to growth from aqueous solution on to 
cleavage surfaces. No x-ray structure data were available at the time, 
but Barker noticed that absence of orientation was always associated 
with large discrepancies between the molecular volumes of the substrate 
and the deposit. The work was subsequently repeated by Royer (1928). 
His results (see table 1) indicated that orientation occurs provided the 
misfit is less than about --14°%. Some of Barker’s results were attributed 
to mixed erystal formation. When orientation occurred it did so not only 
on the (100) cleavage surface, but also on any other surface which was 
used. 
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Sloat and Menzies (1931) examined the effect of using various solvents 
and were able to show that in some cases orientation is obtained with 
misfits larger than 14°. All of their results are included in table 1; the 
effect of different solvents is shown in table 2. It was pointed out that 
reduction in the dielectric constant of the solvent improves the probability 
of orientation, the limiting case being growth from the vapour phase 
(dielectric constant of unity). This was interpreted as showing that the 
solvent interferes with the orientation process and that the higher the 
polarizability of the solvent, the greater is its influence. In these 
experiments, the limiting misfit was extended to +30% (RbI on NaCl). 


Table 2. Effects of Solvent upon Orientation 


(Results of Sloat and Menzies 1931) 


‘ 


Solvent and dielectric constant 
Substrate | Deposit] Water |Furfural Methyl | Ethyl | Acetone} From 
alcohol | alcohol vapour 
81 39 32 25 21 1 
- 
KBr O O O O 
RbCl R O O 
: RbBr R O O O 
Spine ag ea 0 0 O O 
NH,I R R R R R O 
(|. Rb] R R R R O 
KI R R O O O 
KCI NH, I R O O O O O 
; RbI R, O O O O 


¢ . . / . . . . 
R indicates random orientation, O indicates parallel orientation, 


As an extension of these studies, Schulz (1951 c, 1952 a) has examined, 
by electron diffraction, thin alkali halide deposits condensed in vacuo 
from the vapour phase. He finds (table 1) that, for all substrate—deposit 
combinations examined, a parallel-oriented growth occurs, independent of 
the misfit. Misfits varying from —39°, (LiF on KBr) to +90°% (CsI on 
Lif’) are tolerated. Ludemann (1954) has carried out similar experiments 
and confirms the results of Schulz. Elleman and Wilman (1948) find other 
orientations of lead sulphide condensed on to heated rocksalt surfaces. 

Schulz attempted to include an examination of growth from solution 
in his work, by allowing the condensed alkali halide layers to come into 

contact with moist air. Recrystallization occurred when the specimen 
was returned to the vacuum of the diffraction camera. In all cases except 
lithium fluoride substrates, the results were complicated by mixed crystal 
formation. 
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The growth of alkali halides which normally have the caesium chloride 
structure is complicated by polymorphic effects, which are discussed 
in § 5.4. 

3.2. Alkali Halides upon Mica 


The growth of alkali halides upon muscovite mica has been studied 
extensively. The first systematic work was by Royer (1928), who examined 
growth from solution. The most commonly observed orientation has a 
(111) plane of the halide parallel to the (001) mica cleavage surface, with 
the [110] or [110] directions in the halide parallel to the mica [100]. The. 
oriented crystallites are easily visible under the optical microscope, and 
appear as (111) based triangular pyramids. The results are presented in 
table 3. Only salts with a lattice constant of greater than 6-2 4 are oriented. 


Table 3 (a). Alkali Halides with the Rocksalt Structure on Mica 
(Orientation as shown in fig. 3) 


Deposit “0 I II U1 zy: Vv VI M 
Lif 4-02 —45 
NaF 4-62 R oy 
LiCl 5-13 R ZS) 
KF 5:34 R O 57) 
LiBr 5-49 R ess 
AgCl 5:55 R O SoM 
NaCl 5-63 R O O O —23 
RbF 5-64 R = 23 
AvBr 5-77 R O 207 
NaBr 5-96 R Ops}: O —19 
Lil 6-00 R —18 
CsF 6-00 R —18 
KCl 6-28 O O O O —14 
Nal 6-46 O O O —12 
RbCl 6-54 O O —ll1 
KBr 6-59 O O O O —10 
RbBr 6-85 O O — 7 
CsCl 6-94 O — 6 
KI 7-05 O O O O O — 4 
CsBr 7:23 O —1 
NH,I 7:24 O O O —1 
Rbl 7-33 O O O O 0 

O indicates orientation of fig. 3. R indicates random orientation. 


Column. headings : 


I Results of Royer (1928) for growth from solution. 
II Results of Deicha (1946, 1947 a, b, c) for growth from solution. 
III Results of Lisgarten (1954) for growth from solution. 
IV Results of West (1945) for growth from the melt. 
V_ Results of Schulz (1951 b) for growth from the vapour and solution. 
VI Results of Pashley (1952 b) for growth from the vapour. 
M Percentage misfit. 


02 
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Royer considered this in terms of the fit between the (111) face of the 
halide and the potassium ion sites in the mica cleavage surface. The fit 
for the case of potassium bromide is shown in fig. 3, the two parallel networks 


Fig. 3 


mica [olo] 
ker (iz) 


Mica [100] 

Ker [iz0]} 
@ POTASSIUM ION SITES IN (e) POTASSIUM IONS 
MICA (OO!) SURFACE IN KBr (Il!) PLANE 


The fit of the main potassium bromide orientation upon mica. 


Table 3 (6). Alkali Halides with the Caesium Chloride Structure on Mica 
(Orientation as shown in fig. 3) 


Deposit 


Po 


TIC] 3-83 O 45 
NH,Cl 3-87 O ) +6 
NH,Br 4-05 O 411 
CsCl 4-11 R +12 
CsBr 4-29 R O 14 
CsI 4-56 R O ) +25 


O indicates orientation of fig. 3. R indicates random orientation. 


Column headings : 
I Results of Royer (1928) for growth from solution. 
Il Results of Lisgarten (1954) for growth from solution. 
IIt Results of Schulz (1951 ¢) for growth from the vapour. 


each having a hexagonal arrangement. As shown in table 3, the misfit is 
always less than 14°% where orientation occurs. The significance of the fit 
is emphasized by the behaviour of potassium chloride, which is a limiting 
case, A misfit of only —1% occurs if the (111) plane is twisted through 90° 
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(see fig. 4) so as to make the [110] or [110] directions of the potassium 
chloride parallel to the mica [010] direction. This abnormal (11 1) orientation 
was in fact observed, in addition to the above. However, other alkali 


halides which could adopt this orientation with a small misfit failed to 
do so. 
Fig. 4 


A Mica [olo] 
kce [oir] 


MICA [100] 
Kce [277] 
@ POTASSIUM ION SITES IN O POTASSIUM IONS IN 
MICA (OO!) SURFACE Kee (i) PLANE 


The fit of the subsidiary potassium chloride orientation upon mica. 


Table 4. Alkali Halides with the Rocksalt Structure-on Mica 
(Orientation as shown in fig. 4) 


Ao 
Deposit nN I ii Ai M 

NaCl 5-63 N O =) 
NaBr 5-96 N O —6 
KCl 6-28 O 8) O —l 
Nal 6-46 N O +2 
RbCl 6-54 N O +3 
KBr 6-59 N O +4 
RbBr 6-85 N N +8 
KI 7:05 N N +11 
NH,I 7:24 N N +14 
RblL 7-33 N N +16 

O indicates orientation as for fig. 4. N indicates no orientation as for fig. 4. 


Column headings : 


I Results of Royer (1928). 
II ‘vesults of Deicha (1947 b). 
IIL Results of Lisgarten (1954). 


Lisgarten (1954) has re-examined these results of Royer. By using 
very clean conditions, he was able to extend the range of orientation 
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(see table 3). Thus sodium chloride is oriented with a misfit of —23%, 
the crystallites having hexagonal bases as shown in fig. 5. This latter 
case had previously been observed by Deicha (1946, 1947 a, c), who took 
precautions to avoid formation of crystallites in solution, so as to ensure 
that nucleation occurred on the mica surface. In addition, Lisgarten 
observed the abnormal (111) orientation for salts other than potassium 
chloride (see table 4), although this was less frequently observed than the 
normal orientation. 


Fig. 7 


MICA [010] 


mica [loo] 


‘ 
O potassium ion “8° fon] 
@ POTASSIUM ION SITES IN KBr (100) PLANE 
SITES IN MICA ® BROMINE ION SITES IN 
(Oo!) SURFACE KBr (100) PLANE 


The fit of the potassium bromide cube orientation upon mica. 


West (1945) has obtained large oriented sections of certain alkali halides 
by crystallizing them from the melt on a mica cleavage surface. In this 
way, crystals as large as 4 cm diameter by 0-4 em thick have been obtained ; 
the results are included in table 3. The salts are all within the Royer range 
of orientation ; West apparently did not try to grow salts outside this 
range. 

Schulz (1951 b) has deposited alkali halides from the vapour phase in 
vacuo, and finds oriented growth up to misfits of —27%, almost double 
that of Royer (see table 3). He also employed his recrystallization method 
(see § 3.1) for studying growth from solution and found that the results 
were essentially the same as growth from the vapour phase. It has also 
been observed that silver chloride and silver bromide orient from the 
vapour phase on mica with large misfits (Pashley 1952 b, see table 3). 

Deicha (1946, 1947 b, c, 1948, 1949) and Lisgarten (1954) have found 
that certain of the halides give rise to square-based crystallites on the 
mica (see fig. 6). These have (100) planes parallel to the mica surface. 
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Royer (1948 a, b) has criticized Deicha’s results, suggesting that the cubes 
from in solution and then fall on to the mica. Both Deicha and Lisgarten 
interpret their results in terms of nucleation on the mica. Two main 
orientations occur, either the [011] or [001] halide directions being parallel 
to the mica [100] axis. In addition, Lisgarten occasionally finds other 
angular positions. The first of the two orientations is illustrated in fig. 7, 
and the misfits involved are given in table 5. 


Table 5. Misfit Values for Observed Cases of the 
Orientation shown in fig. 7 (Lisgarten 1954) 


°,, misfits °% misfit along mica [010] 
3 when two deposit ion 
cress Along mica Along mica distances are fitted with 

[100] [O10] one substrate distance 
RbBr —6 —46 =f 8 
KBr —10 —48 +4 
RbCl —10 —48 4+ 
Nal —12 —49 +2 
KCl —l4 —5l —] 
NaBr —19 —53 —6 
NaCl —23 —56 —11 

| 


3.3. Alkali Halides wpon Calcite and Sodium Nitrate 
Several workers have studied the growth of alkali halides upon cleavage 
surfaces of calcite and sodium nitrate, which are both rhombohedral 
structures. The crystal axes used in the following description are those 
which define the cleavage rhomb. Heintze (1937) obtained oriented 


Table 6. Misfit Values for Observed Orientations of Alkali Halides 
upon Calcite and Sodium Nitrate 


°% misfit on °% misfit on 
Deposit _ calcite sodium nitrate 
(Royer 1937) (Heintze 1937) 
LiCl —21 
LiBr —15 
NaCl —13 i 
NaBr —7 —8 
KCl —2 —3 
Nal oe 0) 
RbCl oe 
KBr +3 +2 
RbBr +7 
KI +10 +9 
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growth of several halides from solution on to sodium nitrate, and Royer 
(1937) obtained similar results with calcite (table 6). The principle 
orientation observed is illustrated in fig. 8; the cleavage plane of the 
substrates consists of an oblique network of ions which are arranged as 
in a rocksalt cleavage face. There is only one set of parallel lattice rows, 
along the substrate [001] or [010] directions, there being two geometrically 
equivalent orientations. The misfits quoted in table 6 refer to these 
parallel rows. Van der Merwe (1949) also defines an angular misfit between 
the two contact planes. In addition to these orientations, Heintze obtained 
growth with (110) halide planes parallel to the sodium nitrate cleavage 
face with NaCl, KCl, KBr, KI and LiCl. 


Fig. 8 


CoCOs[O10} 
NoBr [o10) 


CALCIUM ION SITES IN © SODIUM ION SITES IN 
CALCITE (100) SURFACE NeBr (100) PLANE 
@ BROMINE ICN SITES IN 

NoBr (100) PLANE 


© CO3 SITES IN CALCITE 
(100) SURFACE 


The fit of sodium bromide on calcite (Royer 1937). 


Schulz (1952 b) has studied thin evaporated layers of a number of 
halides upon both substrates, by means of electron diffraction. Growth 
from solution was examined also. The types of orientation which he 
observed are given in table 7, and table 8 lists the cases for which they 
occurred. Again, the contact planes do not have lattice networks of the 
same form and the misfits quoted in table 8 refer to the fit along the [011] 
and [O11] directions in the substrate surface. The fit of the main 
orientation is shown in fig. 9. A surprising feature of these results is that 
the observed growths from solution do not have the same orientation as 
those found by Heintze and Royer. The reason for this is not clear. 
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The observed misfits are fairly small in most cases. Although there is 
the general tendency for an electrically neutral plane of low indices to 
form parallel to the substrate surface, this condition is not always fulfilled 
(Schulz 1952 b). 


Table 7. Observed Orientations of Alkali Halides upon Calcite 
and Sodium Nitrate (Schulz 1952 b) 


Overgrowth plane Overgrowth axis 
Orientation parallel to substrate parallel to substrate 
(100) cleavage surface [O11] axis 

i (100) [O11] 

Ila (110) [001 | 

IIb (110) [110] 
Ila (111) [110] 
ITI (111) [110] 
IV (741) [113] 


/ 


Table 8. Observed Orientations and Misfits for Alkali Halides 
upon Calcite and Sodium Nitrate (Schulz 1952 b) 


Orientation (see table 7) % misfits 
Growth from Growth from 
Deposit | Structure vapour solution CaCO; NaNO; 
CC  ————e aE ee An ee 
; a Along | Along | Along | Along 

NaCl ‘Ae ii il I I 1 20 2 21 
KCl A I I i it +10] —11 +9} —12 
KBr A if I ul I +16 —7 | +15 —8 
CsCl A i I +22 —2 | +21 —3 
CsBr A U I +27 +2 | +26 +1] 
CsCl B IT IIb +2 |) +16 +1] +21 
CsBr B IIb IIb +6 | +15 +6 | +20 
RbI A I i +29 +4 | +28 +2 
RbI A Mila Ila +4 +2 
KI A I i +24 0 | +23 —1 
KI AN Ila, b | Illa, b 0 —l 
NaF A I U i i) 35 20 35 
NaF A Ila —19 —7 

KF A Ila —6 +7 

LiF A IV IV IV IV ; 

CsF A III —15 


A = Rocksalt type structure. B = Caesium chloride type structure. 
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3.4. Alkali Halides wpon Metals 
Bruck (1936) and Goche and Wilman (1939) found extra spots on 
transmission patterns from oriented layers of silver stripped from rocksalt 
surfaces. Bruck interpreted these as due to oriented rocksalt which had 
crystallized from solution on the silver film ; when rocksalt was condensed 
in vacuo on to heated silver, a different orientation was found (table 9). 


Fig. 9 


caco3(0ol] CoC0s(010} 


(9) POTASSIUM ION SITES 
IN KCL (100) PLANE 


CALCIUM ION SITES IN 


© CALCITE (100) SURFACE 


@ BROMINE ION SITES IN © COs SITES IN CALCITE 
KC£(100) PLANE (100) SURFACE 


The fit of potassium chloride on calcite (Schulz 1952 e). 


Johnson (1950 b, c, d, 1951) has used the optical microscope to study 
the growth from solution of halides upon metal substrates, which were 
in various forms : (1) etched polycrystalline surfaces ; (2) etched massive 
single crystal surfaces ; (3) oriented metal layers formed on other single 
crystal substrates. His results are included in table 9. The misfit for 
many of the results is not quoted because it is likely that the halides were 
growing not on the metal surface but on an oxide surface film (Johnson 
1951). The orientation of the halides on (111) silver surfaces is stated to 
be parallel to the substrate ; hence the data given in table 9. However, 
it is also stated that this corresponds to a misfit of 73° in the case of 
rocksalt. This is consistent neither with the parallel orientation, nor with 
the orientation observed by Schulz (see below) for which the misfit is — 20°. 

Deposits of alkali halides upon silver have been studied in detail by 
Schulz (1952 ¢), who has examined both growth from the vapour phase 
and growth from solution. The substrates were formed by depositing the 
metal on to hot rocksalt to give (100) silver surfaces, and upon hot mica 
to give (111) silver surfaces. The results, given in table 9, show that the 
different methods of growth do not always give the same orientation. 
A wide range of misfits is observed. In all of the experiments listed in 


OE 


Study of Epitaxy in Thin Surface Films 199 


table 9 the topography of the substrate surfaces was not well known ; 
the contact planes have been assumed parallel to the surface for the 
purpose of calculating the misfits. 


3.5. Metals upon Metals 


Most studies of the oriented growth of metals upon metals have been 
made on electrodeposits. x-ray work (Wood 1930, 1935) first indicated 
that epitaxy could occur in such cases. Later, electron diffraction methods 
confirmed this for several substrate-deposit combinations. Cochrane 
(1936) used etched single crystal copper surfaces, and found that some 
metals were oriented on them. Finch and Sun (1936) used beaten foils 
of platinum, palladium and gold as substrates. The results are given 
in table 10. Further data were obtained by Finch, Wilman and Yang 
(1947), who studied deposits on etched copper and iron single crystals. 


Table 10. Electrolytically Deposited Metals 


Type Parallel planes Parallel axes % Retsc. 
Substrate of Deposit mist | “5 ce 
surface Deposit |Substrate| Deposit | Substrate 

Cu B Ni (110) (110) [001] (001 | Sse) 
Cu A Ni (001) (001) 100] [100] = 3) 
Cu B B—Co (110) (110) [001 | F001 | elie (1) 
Cu A B—Co (001) (001) [100] [100 —2 3) 

Cu B Ag (110) (110) (001 | (001 e131 (1).(8) 
Cu A Ag (001) (001) [100] [100 Be vee cca) 
Cu C Ag (T1Ty (111) [110] LLO} | e138 1) 
Cu C Au (111) (111) r110] PL10} | 4 Geis 3) 
Cu G Cr (110) (111) [001] pio |, 2h | a) 

| ey ae Beene +42 ‘ 
Fe B Au (001) (110) [100] oor), fhe |t (3) 
4 Ce 149 
Fe B Ag (001) (110) [100] [001] Bey (3) 
Pt D Cu (110) (110) [001] [001] —8 (2) 
Pt D Ni (110) (110) 001} [001 S100 (2) 
Pt D B—Co (110) (110) [001] [001 —9} (2) 
Pd E Cu (001) (001) 100] [100 —7} (2) 
Pa E Fe (001) (001) [110] [100] 45] (2) 
Au E Fe (001) (001) [110] [100] Sse) 
Au E B—Co | (001) (001) [100] [100] “| ~13 | (2) 
Au E Ni (001) (001) 100] L100} ma 145| (2) 
B-brass EF Cu (211) (211) 111] rO1l] 6 (4) 
2 +5 

B-brass F Cut (1010) (211) [0001] [O11] 0 (4) 


+ Abnormal hexagonal close-packed structure. 


Type of surface : 
A . Etched (001). D_ (110) foil. 


~B_ Etched (110). E (001) foil. 
C Etched (111). F  Electropolished (211). 
References : 
(1) Cochrane (1936). (3) Finch, Wilman and Yang (1947). 
(2) Finch and Sun (1936). (4) Takahashi (1952, 1953). 


A study of copper deposits on an electrolytically polished (211) surface 
of B-brass has been made by Takahashi (1952, 1953). In addition to the 
occurrence of oriented face-centred cubic copper, it is found that a deposit 
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with a hexagonal structure also occurs. This fits with a close-packed 
hexagonal copper lattice and corresponds to a remarkably good fit on 
the B-brass surface, as shown in fig. 10. With all of the oriented electro- 
deposits the misfits are fairly small. Where a bad misfit occurs, it is 
accompanied by a low misfit in the perpendicular direction. 


Fig. 10 
cu [oor] E fort] 
cu [o10] B-BRASS [i777] 
@ ZINC ATOMS IN A-BRASS 
@ ATOM SITES IN COPPER (211) SURFACE 
(1070) PLANE @ copper ATOMS IN £-BRASS 
(211) SURFACE 


The fit of hexagonal copper on a §-brass (211) surface. 


Bruck (1936) deposited silver and gold upon each other, from the vapour 
phase, the substrates being oriented metallic films formed by deposition 
on to rocksalt (see § 3.6). He found that oriented growth occurred when 
the substrate was heated above a certain temperature. 


Table 11. Metals Condensed on to a Silver (111) Surface 
(Newman, results to be published) — 


Orientation 


Metal Structure © _———_____ % 
Thin layer | Thick layer misfit 
Cu A P R —l] 
Au A is R 0 
Ni A R+P R —l4 
Pb A a R+P +21 
Tl ee) BE R+P +19 
Sn c A S+R ia 


Structures : 
A, face-centred cubic. B, hexagonal close-packed. C, face-centred tetragonal. 
\ 


Orientations : 
Hexagonal planes parallel. 
(O11) Sn parallel to (111) Ag with [100] Sn parallel to [110] Ag. 
Random. 
Amorphous. 


b> by 2 by 
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More recently, Newman (results to be published) has studied the growth 
of extremely thin metallic deposits condensed on to atomically flat (111) 
silver surfaces maintained at room temperature. His results are given in 
table 11. The orientations were mostly prominent during the growth of 
the first approximately 104 in average thickness (thin layer) and pre- 
dominantly random orientation occurred with thicker deposits. 


3.6. Metals upon Non-Metals 


The epitaxy of metals upon non-metals has been studied almost entirely 
by electron diffraction. The first evidence of orientation was obtained 
by Lassen (1934). Silver was deposited from the vapour phase in vacuo, 
on to a rocksalt cleavage surface and the resulting metal film was detached 
and examined by the transmission method. A parallel orientation of 
silver was indicated. This was immediately followed by an account by 
Lassen and Bruck (1935 a) of the preparation of good single crystal films 
of silver obtained by deposition on to heated rocksalt. The observed 
orientation corresponds to a misfit of —27°%. This caused Royer (1935) to 
challenge the observations and to point out that a misfit of only +3°% 
would result from the orientation (100) silver parallel to (100) rocksalt, 
with the [011] of silver parallel to the [001] of rocksalt, which would 
therefore be expected to occur. However, Lassen and Bruck (1935 b) 
confirmed that their original claim was correct and that the latter 
orientation is not observed, despite the small misfit. This was confirmed 
by Kirchner and Lassen (1935), who used the reflection technique. The 
same result has since been obtained by many workers including Finch 
and Wilman (1937), Kirchner and Cramer (1938), Goche and Wilman 
(1939), Uyeda (1942), Shirai (1943 b), Johnson (1950 a) and Goswami 
(1954). Bru and Gharpurey (1951 a,b) and Raether (1951) also find 
parallel orientation of silver on (110) and (111) surfaces, where misfits 
are also —27%. 

Bruck (1936) has studied the growth of a number of metals on rocksalt 
and finds a marked dependence upon the temperature of the substrate. 
Above a certain critical temperature (the epitaxial temperature) good 
orientation was obtained ; below this temperature, the growth included 
at least some randomly oriented deposit. Bruck’s results for the face- 
centred cubic metals are given in table 12. A characteristic of many of 
the above experiments is the occurrence of pronounced twinning on the 
(111) planes. This is discussed further in §§ 5.1 and 6.1. 

The results for body-centred cubic metals are given in table 13. This 
includes data by Shirai, who has carried out many experiments by means 
of the transmission technique. In several cases the substrates were 
preheated to a temperature of several hundred degrees Centigrade before 
the deposition was carried out at some lower temperature. This appeared 
to influence the results. In addition to the orientations quoted by Bruck 
(1936), Shirai found that orientations with (111), (210) and (331) planes 
parallel to the surface occurred, the results being erratically dependent 


Table 12. Observed Cases of Parallel Orientation of 
Face-Centred Cubic Metals upon Alkali Halides 


: Epitaxial A References 
Substrate Deposit So ie ight 
Ni 370 —38 (1), (4), (2) 
B—C, +540 —37 (1) 
Cu 300 —36 (1) 
NaCl Pdt 250 —31 (1), (9) 
Alt 440 —28 (1) 
Aut 400 —28 (1); (2), (6) 
Ag 150 —28 (1), (3) 
Ni 410 —44 (1) 
Kn Pd 250 —38 (1) 
BC Ap 150 35 (1), (3) 
L Pt —38 sa2) 
KBr Ag 115. —38 (3) 
KI Ag 80 —42 (3) 
NaBr Ag 80 —32 (3) 
References : 
(1) Bruck (1936). (5) Fordham and Khalsa (1939). 
(2) Shirai (1948 a). (6) Trillat, Terao, Tertian and 
(3) Shirai reve b). Gervais (1955). 
(4) Collins and Heavens (1954). (7) Thirsk (1950). 


+ (100) Pd parallel to (100) NaCl, with [011] Pd parallel to [001] NaCl, 


observed also (ref. (5)). A! 
t (111) metal parallel to (100) NaCl, with [110] metal parallel to [011] NaCl, 


observed at lower temperatures (refs. (1) and (2)). 


Table 13. Observed Cases of Orientation of Body-Centred Cubic Metals 
upon Alkali Halides 


°% misfit 
Substrate Deposit | ————————-———————_|_ Reference 
Orientation | Orientation 


References : 
(1) Bruck (1936). (3) Shirai (1939). (5) Shirai (1941 b). 
(2) Shirai (1938). (4) Shirai (1941 a). (6) Shirai (1937). 
Orientations : 


[ (100) metal || (100) halide, with [011] metal || [001] halide. 
II (110) metal || (100) halide, with [001] metal || [001] halide. 


Table 14. Metals upon Mica 


Epitaxial Plane Axis 
Metal temperature | parallel to | parallel to oe References 
(76) surface mica [100] misfit 
( 150 (111) [110] —44 (1), (2) 
"oe (111) [211] ae (1), (2) 
250 (100) [001] ae (2) 

. 450 (111) [110] —44 (1) 
a (111) [211] aA (1) 
- (111) [110] ej aa) 

470 a aT 
: (111) [217] ot a) 
References : 
(1) Rudiger (1937). (2) Thirsk (1950). 


Table 15. Metals upon Calcite (Rudiger 1937) 


‘ Epitaxial Plane Axis 
Metal temperature | parallel to parallel to i 
(6) calcite (100) | calcite [001] misfit 
surface 
470 (111) [110] =i 0) 
Ag 470 (111) (21T] 22 
510 (100) [011] 10 
360 (111) [110] 10 
4 360 (111) [211] —22 
aa 470 (100) [011] —10 
510 (110) [112] —22 
400 (111) [110] —15 
Pd 400 (111) [211] —26 
420 (100) (011) a5 
L 490 (110) (112 —26 


Table 16. Metals upon Fluorspar (Rudiger 1937) 


Epitaxial % misfit 
Metal temperature © |_——$——____—__ 
(°c) Orientation I | Orientation II 

ae 2) eee aera 

Ag +500 25 ae 

: = 30 

Au 380 —25 su 

. +23 

Pd 400 —29 "59 


Orientations :_ = 
T (111) metal || (111) CaF, with [110] metal || [110] CaF, 
II (111) metal || (111) CaF, with [211] metal || [110] CaF, 


P.M. SUPPL.—APRIL 1956 
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upon temperature. This behaviour was found for chromium (Shirai 1939), 
molybdenum (Shirai 1941 a, b) and iron (Shirai 1937, 1938), which was 
found to behave similarly upon KCl, KBr and KI. 

Oriented growth of metals upon mica, calcite and fluorspar te also 
been found (see tables 14, 15, 16). These are characterized, in general, by 
the occurrence of mixtures of orientations which are temperature 
dependent. The misfits are often quite high, but can be considered as 
smaller if two deposit atom distances are compared with one corre- 
sponding substrate distance. The misfits quoted in table 15 refer to 
the calcite [001] direction only, since no other parallel lattice rows occur. 
In addition to the results quoted in table 14, Votava (1953) has found 
oriented growth on mica of platinum formed by a method due to Pigeon 
(1890). The metal crystallites in the form of hexagonal plates were 
visible under the optical microscope. 


Table 17. Metals upon Mineral Sulphides 


i 
Temperature | Parallel planes Parallel axes % 
Substrate | Deposit | conditions |—————~—————|_————_——| misfit | References 
(°c) Sub. Dep. Sub. Dep. 
Pt > 250 (0001) | (111) | [100] | [110] | +4 (6) 
Ag > 20 (0001) | (111) | [10I0] | [110] +9 | (1), (2). (6) 
Mos, Cu > 50 (0001) | (111) | [1010] | [110] —3 (6) 
Ni > 120 (0001) | (111) | [1010] | [110} —6 (6) 
Au > 80 (0001) | (111) | [1010] | [110] +9 (6) 
Pbs { Cu 100 to 250 (100) | (110) (011] | [111] i (4) 
Ag > 100 (100) | (100) [001] | [001] | —31 (1), (5) 
Au > 120 (110) | (110) (001] | [001] | —25 (3) 
Ag > 80 (110) | (110) {001} | [001] | —25 (1), (2) 
Ni > 120 QUAD le i [001] | [110] ist (3) 
rte Cu <190 (110) | (310) | [1T0} | [001] Ave (3) 
Cu 190 to-350 (110) | (110) (001] | [001] | —33 (3) 
‘ 850) to 5 7 —53 (3) 
Cu 350 to 500 (110) (111) [001] [110] +15 
Fes, Ag > 150 (100) | (100) | [oor] | poor: | —24 (1) 
References : 
(1) Uyeda (1940). (3) Miyake and Kubo (1947 a). (5) Kubo and Miyake (1948). 
(2) Uyeda (1942). (4) Miyake and Kubo (1947 b). (6) Kainuma (1951). 


The deposition of metals on to mineral sulphide surfaces obtained by 
cleavage has been studied by a number of workers. Uyeda (1940, 1942) 
has studied silver on galena, molybdenite and zinc blende, and Miyake 
and Kubo (1947 a, b) have studied deposits upon galena. In the latter 
experiments the results were complicated by effects attributed to the 
diffusion of the metal into the galena, and to the occurrence of chemical 
reactions which lead to the formation of oxides and sulphides of the 
deposit metal, as well as of metallic lead. The orientation of these reaction’ 
products is considered in § 3.7. Where oriented metallic deposits were 
formed, the results are given in table 17. Kubo and Miyake (1948) have 
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also studied the growth of metals upon zinc blende. No misfit is quoted 
for the orientation of gold given in table 17, because it is stated that the 
(110) gold planes were not exactly parallel to the surface. For the growth 
of metals upon molybdenite (Kainuma 1951), a (111) plane of the metal 
is parallel to the hexagonal basal plane of the substrate ; mostly there is 
double positioning, but in the case of nickel and gold, single positioning 
occurs under some conditions. 


3.7. Chemically Grown Layers 


With the exception of a few x-ray studies, all of the results for the 
occurrence of preferred orientation in chemically grown layers have been 
obtained by electron diffraction. These results are given in table 18. 
Several types of reactions are included. In some cases the contact planes 
between the substrate and the overgrowth are not known, so that a 
significant value of the misfit cannot be determined. In other cases, the 
contact planes have been determined by special techniques (see § 2.3.3), 
or are fairly certain as a result of the type of substrate used. For these, 
a significant misfit value can be calculated, and has been included in 
table 18. These values vary between wide limits. 


Fig. 12 


438r[00}] Ag [ITO] 


AgBr [110] Ag [Oo] 
C:) SILVER IONS IN AgBr (110) PLANE 


@ ATOM SITES IN SILVER (110) SURFACE 
e BROMINE IONS IN AgBr (110) PLANE 


The fit of silver bromide on a silver (110) surface. 


The chemically grown deposits have provided several clear examples: 
of the dependence of the orientation upon the contact planes. Thus for 
silver halides upon silver, the orientation on the (100), (110) and (111) 
faces are quite different. It was originally found that silver bromide on 
the (111) face of a massive silver crystal is mainly in a fibrous orientation 
(Pashley 1952 a). However, more recently the author has carried out 
many experiments with substrates consisting of atomically flat highly 
oriented silver layers on mica, which have (111) planes parallel to the 
surface. When fresh surfaces are used, the orientation is not fibrous, but 
is very well defined, with the (111) planes of the bromide parallel to the 
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surface and [110] or [110] bromide directions parallel to the silver [110] 
axis (double positioning). In some cases preferred single positioning 
occurs ([110] AgBr only parallel to Ag [110]). This is illustrated in fig. 11 ; 
double positioning would give the same array of spots, but the intensity 
distribution in the pattern would be symmetrical about the central line. 
These orientations correspond to a +41°%, misfit between two parallel 
hexagonal networks. The exact nature of the orientation appears to 
depend upon thickness and is still being investigated. In contrast, the 
orientations on the (100) and (110) silver surfaces correspond to a zero 
misfit (see fig. 12). 

An interesting effect is found for the growth of zinc oxide on a heated 
zinc blende surface. The hexagonal zinc oxide always grows with its 
[0001] axis nearly parallel to a [111] type axis of the sulphide, the exact 
orientation being slightly adjusted so that a particular lattice plane is 
exactly parallel to the surface, which has in all cases been parallel to a 
closely packed plane of the sulphide lattice. Since there are four distinct 
[111] type axes in the zine blende, there are four possible orientations 
on any surface. Not all of these possibilities occur on a given face, the 
results depending upon the direction of the piezo-electric axis (Uyeda, 
Takagi and Hagihara 1940, 1941). 

An effect which appears to be unique in studies of growth on single 
crystal substrates has been observed by Frisby (1948). Various electro- 
polished copper surfaces which contain a [110] axis, but which are not 
necessarily parallel to any lattice plane, have been oxidized in water 
at 100°c. It is found that a (110) plane of the cuprous oxide always grows 
parallel to the surface, with the [110] direction parallel to the copper 
[110] axis. This appears to be the only known case, excepting ‘ fibrous ’ 
orientations, where the arrangement cannot be defined in terms of two 
parallel lattice networks which have at least one pair of parallel lattice 
rows. 
Most of the x-ray work has been carried out by Schwab (1942, 1947 a, b,c). 
Because of the uncertainty of the determination of the contact planes in 
this work, the results are not included in table 18. Oriented growths have, 
however, been established for a number of reactions, including halides 
upon silver, and silver and thallium halides upon each other by ‘exchange’ 
reactions. 

According to the table published by van der Merwe (1949), magnesium 
oxide is oriented upon magnesium. The reference quoted is that of Finch 
and Quarrell (1933). This is an incorrect interpretation of the observed 
facts, and only followed from an unconfirmed remark concerning pseudo- 
morphic growth (see § 4.2). There appears to be no valid evidence that 
the magnesium oxide is oriented upon magnesium. 


3.8. Orientation during Decomposition 


In addition to the occurrence of epitaxy in thin chemically grown layers, 
it is also found that orientation can occur during decomposition processes. 
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Berry (1949) and Berry and Griffith (1950) irradiated single crystals of 
silver bromide with ultra-violet. x-ray diffraction examination showed 
that silver in parallel orientation (two parallel cubic lattices) grows well 
inside the silver bromide, while electron diffraction examination showed 
that other orientations occur at the surface. Pashley (1951 b, 1952 ¢c) 
has decomposed thin oriented layers of silver halide by electron bombard- 
ment in the electron diffraction camera (reflection technique). He finds 
that the silver is in a well-defined parallel orientation when the chloride 
and bromide are used, but that less marked orientation occurs with the 
iodide. Prominent (111) twinning is also observed. Although Trillat 
(1951, 1952) failed to find such orientation when thin oriented silver 
bromide films were examined by the transmission technique, Pashley 
(1952 d) has obtained orientation with such specimens of both chloride 
and bromide. The pattern of fig. 13 illustrates this. The square array 
of very intense spots is due to unchanged silver chloride, and the other 
weaker spots arise from silver in parallel orientation, including extra 
spots due to secondary diffraction. The absence of orientation in Trillat’s 
films was explained as due to a gradual disorientation of the bromide, 
before decomposition, resulting from heating effects. 

The decomposition of small crystals of silver azide has been studied by 
Sawkill (1955), who finds that much silver is formed with its cube planes 
parallel to the (100) plane of the azide, which has a body-centred ortho- 
rhombic structure. A [110] silver direction is parallel to the azide [001] 
direction. . 


§ 4. Spactna MEASUREMENTS AND ‘ BASAL PLANE PSEUDOMORPHISM ” 
4.1. The Early Work and its Original Interpretation 


The term ‘basal plane pseudomorphism’ was introduced by Finch and 
Quarrell (1933) in the explanation of certain of their results concerned 
with sputtered metal films. It was found that zine deposited on to a 
quartz substrate rapidly became covered with an oxide layer. Measure 
ment of the electron diffraction rings obtained by the reflection technique 
revealed that a certain number of the rings fitted with neither metallic 
zinc nor zinc oxide. These extra rings were interpreted as arising from 
a new structure of the oxide, which is of the normal ZnO type, but which 
has a modified cell shape. The a-axis is reduced to equal that of the 
hexagonal metallic zine lattice, while the c-axis is increased so as to 
maintain an approximately unchanged cell volume. The occurrence of 
the new cell was attributed to the influence of the metallic zine structure. 
It was considered that the oxide grew in contact with the basal (0001) 
planes of the zine and in parallel orientation, such that its basal plane 
was constrained to fit that of the zinc. This * pseudomorphic ’ structure 
was believed to grow to a considerable thickness before recrystallization 
produced the normal structure. Later, Finch and Quarrell (1934) studied 
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by the transmission technique thin skins of zinc oxide which had been 
lifted on a gauze from the surface of molten zinc exposed to air. It was 
claimed that the resulting ring patterns also indise red the presence of this 
pseudomorphic structure. 

Other similar cases have been reported. Finch and Quarrell (1933) 
deposited a layer of aluminium on to a previously sputtered layer of 
platinum, and concluded that the aluminium had formed in a new 
tetragonal structure. The basal plane of the new cell is supposed to fit 
the cube plane of the platinum. A layer of magnesium oxide (Finch and 
Quarrell 1933) formed by oxidation of sputtered magnesium was found 
to have abnormal spacings which were also attributed to pseudomorphic 
effects. 

‘Cochrane (1936) found evidence for modified spacings in electrodeposited 
layers of nickel and cobalt on single crystals of copper. With a fast rate 
of deposition he obtained randomly oriented films of nickel, which gave 
rise to a pattern of rings which passed through the diffraction spots 
attributed to the underlying copper. It was concluded that all the spacings 
of the nickel lattice were abnormal and equal to that of the copper. With 
slow rates of deposition, Cochrane found that parallel orientation of nickel 
and cobalt occurred, and that this was accompanied by repeated twinning 
on the (111) planes, which gave rise to extra features on the diffraction 
patterns. Measurements indicated that the nickel and cobalt had spacings 
equal to that of the copper. Although there was no change in cell shape 
in these cases, it did appear to suggest some form of influence of the 
substrate on the lattice spacings of the overgrowth. . 

Miyake (1938) obtained anomalies in the spacing of Sb,O, formed on 
stibnite (Sb,S;) which he has attributed to the influence of He stibnite 
substrate. 

Clark, Pish and Weeg (1944) have studied the electroplating of several 
metals on to various metal sheets. In some cases, the absence of x-ray 
diffraction patterns from the thinnest layers is interpreted as due to a 
three-dimensional pseudomorphic effect. 


4.2. The Validity of the Early Results 


Despite the above evidence for the occurrence of basal plane pseudo- 
morphism, many studies of even thinner layers of well-oriented deposits 
have failed to produce further evidence. Since the concept of ‘ basal 
plane pseudomorphism’ is frequently put forward in textbooks and 
_ reviews, it is of importance that the position should be thoroughly 

examined, in an attempt to discover why more evidence of the effect 
is not obtained. It is natural that this investigation should start with 
an examination of the validity of the earlier results. 

With the exception of the work of Cochrane (1936) and Miyake (1938), 
none of the early experiments was carried out on single crystal substrates. 
Thus, in the case of the oxide on zine, it was not even established that the 
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oxide was oriented, although this was an implicit assumption. In addition, 
only growth on hexagonal (0001) planes of the zinc was considered, 
although other types of facets, upon which the occurrence of different 
pseudomorphic structures would be expected, would certainly be present 
on a polycrystalline surface. More recently, therefore, other workers have 
investigated this aspect. Raether (1950), Lucas (1951, 1952), Ehlers and 
Raether (1952) and Ehlers (1953) have studied the growth of zine oxide 
on single crystals of zinc. The first two papers deal with the growth which 
occurs upon a freshly cleaved surface, as soon as it is exposed to air. The 
oxide forms almost instantaneously, so that examination of the surface 
by electron diffraction gives rise to oxide diffraction patterns, but no 
patterns from the underlying metal. It is always found that the zine 
oxide is oriented parallel to the zinc, and that it has its normal structure 
and lattice spacings. No evidence was found for the pseudomorphic 
structure. Attempts have been made to study the initial stages of 
oxidation, but the extreme affinity between zine and oxygen has so far 
limited this aspect. Raether estimates that his oxide film is at least 
20 to 40 A in thickness. By cleaving a specimen and immediately plunging 
it into benzene, Lucas (1951) was able to reduce the thickness of the oxide 
sufficiently to allow the zinc substrate pattern to be obtained. The 
specimen was mounted with its surface still covered by benzene, which 
evaporated as the diffraction camera was evacuated. The film thickness 
under these conditions was estimated to be certainly less than 25 A, it 
being difficult to make any accurate determination. Again the oxide had 
its normal structure. Although thinner layers could probably be studied 
if controlled growth were carried out on a specimen cleaved ina vacuum 
the main conclusion follows quite clearly without such experiments. The 
most favourable condition for the growth of the pseudomorphic form 
reported by Finch and Quarrell would undoubtedly be on the cleavage 
(0001) face of the metal. The experiments of Lucas and Raether show 
that such a structure does not occur, except possibly during the growth 
of the first two or three atomic layers only. But if the structure did occur 
in the first two or three layers only, it would not have been observed by 
Finch and Quarrell, since their experimental conditions would not have 
allowed the detection of such thin layers as readily as in the experiments 
of Raether (1950) and Lucas (1951). It therefore follows that the extra 
rings observed by Finch and Quarrell were not due to a pseudomorphic 
form of zine oxide occurring as a result of epitaxy. Thus the existence 
of such a pseudomorphic growth is reduced to that of theoretical 
speculation and has no valid experimental foundation. 

Experimental checks of the case of aluminium on platinum are not so 
easy. Although subsequent reports of this experiment have described it 
as growth on to ‘a erystalline platinum substrate in (100) orientation ’ 
(Quarrell 1941), this is an unjustified description. The actual substrate 
consisted of a layer of platinum sputtered on to a quartz surface. This 
layer was randomly oriented, A layer of aluminium was then evaporated 


Study of Epitaxy in Thin Surface Films 213 


on to the platinum from a filament, the heating from which caused 
recrystallization of the platinum during the deposition, to give an oriented 
film (see Finch and Quarrell 1933, p. 413). A repeat of this experiment 
under more controlled conditions has probably not been carried out 
because of the difficulty of preparation of a suitable substrate. However, 
a close examination of the evidence of Finch and Quarrell (1933) reveals 
a completely different interpretation of the observations. 

Finch and Quarrell appear to have treated their a pattern as one of 
continuous rings and to have used only the radii of the rings for the 
assignment of indices. They do not appear to have correlated the indexing 
with the distribution of the ares of the pattern (see their fig. 9). However 
their published pattern is a contact print of the actual photograph and 
sufficiently accurate measurements can be made on this print to allow 
the unambiguous identification of the original measurements for almost 
every arc. In this way, the assignment of indices by Finch and Quarrell 
can be studied. This indexing is found to be inconsistent, as has already 
been pointed out by Shishakov (1952). 

During the recrystallization process, the platinum substrate film orients 
in a fibrous manner, with its (100) planes parallel to the surface (see Finch 
and Quarrell 1933, fig. 12). There are two possible ways for the growth 
of the tetragonal aluminium. 

Firstly, the (001) aluminium planes are parallel to the platinum surface 
(100) planes. This would give an (001) aluminium fibrous pattern, which 
would be a slightly distorted form of the published platinum pattern 
(see Finch and Quarrell 1933, fig. 12). This is clearly not the case, since 
the pattern obtained after the deposition of the aluminium shows pairs 
of ares in place of the single arcs given by the platinum. The second 
possibility is that the aluminium (001) planes are parallel to either the 
(100); (010) or (001) platinum planes. This would correspond to the 
aluminium growing in a ‘ pseudomorphic ’ manner on all of the cube-type 
‘ platinum facets. This type of growth would give rise to the pattern 
plotted in fig. 14. It will be seen that the pattern is of the same basic 
form as a simple (001) fibrous pattern, but that many of the spots now 
become split into pairs or triplets. Part of this splitting would be obscured 
by the arcing of the patterns, but some of it would still remain observable. 
Comparison with Finch and Quarrell’s pattern shows certain differences 
in form and many differences in indexing. Although the absence of 
triplets (e.g. the (420), (402) and (204) set) can possibly be explained as 
due to insufficient resolution (e.g. (402) and (420) unresolved), the 
observation of a doublet where a single arc is predicted (e.g. the (222) 
set) cannot be explained in this way. Also, the variation in the spot 
separations from pair to pair does not agree. Thus, in accord with the 
conclusions of Shishakov (1952), it is necessary to reject Finch and 
Quarrell’s interpretation of the pattern. 

The pattern may be completely and satisfactorily interpreted as due 
to a mixture of platinum and aluminium with their normal cubic structure 


214 D. W. Pashley on the 


and lattice spacings. This interpretation is illustrated in fig. 15. The 
form of the pattern appears to be identical with the observations. The 
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The indexed pattern which would be obtained from Finch and Quarrell’s 
(1933) supposed pseudomorphic tetragonal aluminium layer. A super- 
position of (100) and (001) fibrous patterns. 
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The indexed pattern which would be obtained from a mixture of normal 
aluminium and platinum in (OOL) fibrous orientation. 


agreement with measurements is given in table 19. Finch and Quarrell’s 
measurements have been used, but their calculated spacing values have 
been changed by + 2-8%, corresponding to a wavelength error of +2:8%, 
to give the best agreement with the normal aluminium and platinum 
spacings. This adjustment is justified because Finch and Quarrell claim 
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an absolute accuracy of only about 2°% for their measurements of the 
wavelength. This reinterpretation agrees with that put forward by 
Shishakov (1952), except for a few minor differences in indexing. However, 
it is considered that Shishakov’s further deduction that the platinum has 
a slightly abnormal spacing is unjustified in view of the limited accuracy 
of the measurements. 


Table 19. Interpretation of Pattern from Aluminium—Platinum Layer 


Data from Finch and Quarrell (1933) Reinterpretation of the results 
(Akl) Aluminium Platinum 
r tetragonal dA da 
(em) aluminium | (A=0-0564) | (A=0-0580) | (hkl) dA (hkl) da 
1:17 (111) 2°29 2°35 (CL )) 2-33 (111) 2-267 
1-33 (002) 2-01 2°07 (200) 2-02 (200) 1-96F 
1-95 (220) 1°37 1-41 (220) 1-42 (220) 1387 
2-25 (113) 1-19 1-22 (311) | 1-22 
2-31 (131) 1-16 1-19 (S11) oh-3-18 
2-36 (222) 1-14 1-16 (222) | 1-16 
9-41 1-14 (222) 1-13 
2-75 (400) 0:98 1-00 (400) 1-01 (400) 0-98+ 
2-98 (133) 0-90 0-92 (331) 0-93 
3:05 (204) 0-88 0-90 (420) 0-90 (331) 0-90 
3-12 (240) 0:86 0:88 (420) 0:88 
. 3°34 (224) 0-80 0-82 (422) 0-82 
3:4] (242) 0-79 0°81 (422) 0-80 
3-54 { ie 0-76 0-78 (511) | 0-78 
3-61 (151) 0-74 0-76 (511) | 0-76 
3°90 (440) 0:69 0-71 (440) 0-71 
397 (135) 0-67 0-69 (440) 0-69 
(153) é : 
4-08 (244) 0-66 0-68 (531) 0-68 
(351) A) 
4-15 (600) 0-65 0:66 (531) 0-66 
(442) 
4-35 (260) 0-62 0-63 (620) 0-64 
sa) { ee 0-61 0-62 (620) | 0-62 
4-50 (353) 0-60 0-61 (622) 0-61 5 
4:60 (262) 0:58 0-60 (622) 0:59 
4:89 (406) 0-55 - 0-56 (711) 0:56 a 
5-00 (460) 0-54 0-55 (711) 0:55 
» 


tAluminium and platinum ares unresolved. 


This interpretation of the pattern is perfectly reasonable. The aluminium 
is deposited on to randomly oriented platinum, which recrystallizes during 
the deposition because of heating effects. A mixture of aluminium and 
platinum is to be expected. It is therefore concluded that the evidence 
for a ‘ pseudomorphic ’ tetragonal form of aluminium is false. 

The experiments of Cochrane (1936) concerning nickel and cobalt 
deposits have been repeated in this laboratory by Newman (1956). 
He finds no change in spacing of nickel deposited at a fast rate 
(see § 4.1). The pattern interpreted by Cochrane as due to nickel 
rings passing through copper spots was, in fact, a pattern of nickel rings 
passing through nickel spots. Apparently some of the nickel was oriented, 
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even at the faster rates of deposition, giving rise to a nickel spot pattern, 
and sufficient nickel had been deposited to obscure the copper pattern. 

For the slow rates of deposition, Newman confirmed Cochrane’s result 
that the spacing of the deposited layer was that due to copper and not that 
of either nickel or cobalt. The electrolytic bath was being run under 
conditions for which no deposition would be expected, because the cell 
voltage was not above the overpotential. However, both Cochrane 
(private communication) and Newman obtained blackish-looking deposits. 
This deposition occurred because of unclean anodes and it was found by 
Newman that no deposition occurred when a carefully cleaned anode was 
used. This was checked, in the case of the cobalt bath, by the use of a 
radioactive cobaltous sulphate solution. It was found, by a tracer 
technique, that less than 1A of cobalt was deposited, no matter how 
long the cell was left running, for the low current densities. However, 
an interesting result followed. The specimens treated in this way often 
showed considerably changed diffraction patterns, although no cobalt had 
been deposited. The patterns were very similar to those obtained from 
the specimens which showed the blackish deposits, prominent twinning 
spots and streaks being observed. This has been interpreted as resulting 
from the dissolution and subsequent redeposition of copper on to the 
specimen. This evidence gives a very plausible explanation of the spacing 
measurements on the deposited layers (Newman, to be published). The 
cobalt deposition occurs at a very slow rate, and is accompanied by the 
redeposition of copper. This latter would tend to deposit preferentially 
on the specks of deposited nickel or cobalt. This would result in the 
deposited nickel or cobalt being covered by a layer of copper, which would 
prevent its being observed by electron diffraction. It is also possible that 
small regions of nickel-copper or cobalt—copper alloy are formed. Even 
if these alloy regions did contribute to the diffraction pattern, their 
spacings would be indistinguishable from that of copper, unless they 
contained more than about 20%, of the nickel or cobalt. Such an 
explanation of the observed spacing measurements seems far more 
plausible than that of the formation of metallic nickel or cobalt with 
lattice spacings equal to that of copper. 

The other remaining supposed evidence of pseudomorphic growth (see 
§ 4.1) is even less convincing. Miyake’s (1938) identification of an Sb,O, 
layer on Sb,Sz is complicated by the possibilities of other oxide formation. 
Finch and Quarrell’s (1933) magnesium oxide layer gave a pattern of a 
few rather broad rings. If a wavelength error is allowed, most of these 
rings may be fitted to cubic magnesium oxide. The weakness of the 
evidence of Clark, Pish and Weeg (1944) is that their evidence is of a 
negative character. Certain electrodeposited layers of about 500.4 mean 
thickness gave no X-ray patterns, whereas 1000 A layers did. There are 
several other possible interpretations of these results, besides the ‘ pseudo- 
morphic ’one. In any case, it is difficult to visualize (to quote one example) 
how a pseudomorphic zine structure could be identical, as regards all 
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spacings, with that of steel. In view of the lack of detailed evidence in 
these three cases, they are rejected in that they do not supply conclusive 
evidence in favour of the hypothesis that ‘ pseudomorphic ’ growths occur. 
From the above examination of all of these early results it follows that 
m no case is the claim of the existence of a pseudomorphic structure 
justified. None of these experiments supports the model of pseudomorphic 
growth in epitaxy. 


4.3. The Extent to which Pseudomorphism can be Investigated 
Haperimentally 


Since the early evidence on pseudomorphism is invalid, it is necessary to 
consider to what extent the phenomenon can be investigated experimentally. 
Many investigations of oriented layers which are less than 10 4 in average 
thickness have been carried out in this laboratory and elsewhere (see, for 
example, Schulz 1951 b, 1952 a, Uyeda 1942). In none of these experiments 
has positive evidence of pseudomorphism been found. 

Electron diffraction can be used to analyse single atomic layers on flat 
surfaces, if the growth occurs in the form of monolayers (see § 2.3.6). 

‘The experiments which most nearly fulfil these conditions are those of 
Schulz (1951 b, 1952 a) concerning the growth of alkali halides upon 
other alkali halides and upon mica. When the misfit between overgrowth 
and substrate is less than about 10°, a deposit of less than one atomic 
layer in average thickness does not form independent widely separated 
nuclei, but collects into patches of some unknown size, but not more than 
an atom or two in thickness. Although Schulz (1951 b) originally 
suggested that the ‘monolayer’ patches had spacings equal to that 
of the substrate, he appears to have modified his views, and has since 
stated (Gomer and Smith 1953, p. 375): “...There was no evidence to 
show that the natural bulk spacing of the deposit changed to match that 
of the substrate ; the evidence (what there was) would suggest that the 
deposit always retained its natural bulk spacing...”. In addition, 
Blackman (1952) has pointed out that (111) alkali halide monolayers (on 
mica) would be of a special type, whose natural spacing would not be that 
of the bulk material. The detailed interpretation of the almost monolayer 
patterns is made difficult by their diffuse nature, which limits the accuracy 
of measurement to a few per cent. Further ambiguities could also arise 
as a result of secondary diffraction, which could shift the centre of a broad 
diffuse streak. This latter effect is illustrated in fig. 16. In addition to 
the primary diffuse streaks from the deposit (assumed, for this argument, 
to have its natural bulk spacing), secondary streaks could arise as a result 
of the substrate diffractions acting as new sources. The spots along the 
central line would give secondary patterns coincident with the primary 
patterns. The spots along the rows A of the substrate pattern, for example, 
would, if they acted as secondary sources, give rise to the extra diffuse 
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streaks indicated. These would merge along the lines B with the primary 
streaks to give broad diffuse streaks which would appear to have the same 
spacing as those of the substrate. In this way, a deposit which is a few 
per cent different in spacing from that of the substrate could give a diffuse 
streak pattern which would erroneously suggest a conformity in spacing. 
We have made some attempts in this laboratory to investigate in more 
detail the monolayer region found by Schulz, by using the continuous 
examination technique described in § 2.3.5. However, no further evidence 
has been obtained. With mica substrates, the substrate pattern is not 
sufficiently well defined to allow a conclusive study of the cases of small 
misfit (10°, and less), because of charging effects with the mica. With 
alkali halide substrates it has been found that the continuous examination 
technique is not very successful because the electron beam affects the 
substrate in some manner which seriously influences the growth of the 
deposit. Our general experience so far would support the conclusion of 
Schulz which is quoted above. There is certainly no doubt that, once the 
deposit becomes more than about one atomic layer in average thickness 
separated crystallites are formed and that these have their natural bulk 
spacing. 
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When the deposit grows in such a manner that it has formed widely 
separated crystallites when it first gives an observable pattern, the natural 
bulk spacing is always observed. At this stage, most of the crystallites 
will already have grown to a certain thickness, which will vary con- 
siderably from case to case. The value of the deposit lattice spacings at 
any earlier stage of their growth cannot be measured. 
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4.4, Conclusions 


From the evidence available to date, it seems reasonable to conclude 
that pseduomorphism, in the sense in which it was first postulated by 
Finch and Quarrell, does not occur. Whether it occurs in a much more 


restricted sense is open to speculation, but certainly no positive evidence 
of this has been found. 


§ 5. EXPERIMENTAL EVIDENCE CONCERNING THE MopE oF GRowTH 
OF AN ORIENTED LayER 


Before any comprehensive theory of epitaxy can be formulated, a 
certain amount of knowledge of the mode of growth of an oriented layer 
is essential. Although optical microscope studies can give useful infor- 
mation, most of the existing direct evidence concerning the earliest stages 
of growth has resulted from electron diffraction studies. 


5.1. Monolayers and Nucleated Deposits 


One of the most important questions is whether the initial deposit 
grows in the form of a two-dimensional monolayer, or whether three- 
dimensional nuclei or crystallites are formed initially. It is possible to 
obtain the most systematic evidence in the case of growth in vacuo from 
the vapour phase, because this may be carried out inside a diffraction 
‘camera. For reliable results, it is necessary that an atomically flat substrate 
surface be used. Schulz (1951 b, 1952 a) has studied the growth of alkali 
halides upon other alkali halides and upon mica. He finds that the initial 
deposits are in the form of three-dimensional crystallites when the atomic 
misfit is greater than a certain value (20% for alkali halide substrates 
and 10°, for mica substrates), but that the growth of ‘ monolayers ’ 
occurs for smaller misfits. The evidence for these “ monolayers ’ derives 
from the observation of streak patterns from deposits which are less than ~ 
one molecular layer in average thickness. Because of the weakness and 
diffuse nature of these streaks, detailed interpretation is difficult. Schulz 
deduces that the deposits occur in patches of not more than an atom or 
two in thickness, but of unknown lateral extent (see Gomer and Smith 
1953, p. 375). It is certainly quite clear that the deposit material spreads 
over the substrate surface far more in these cases than it does for cases 
of larger misfit. This kind of effect was observed as long ago as 1907 by 
Barker (1907, 1908), who obtained so-called zonal growths (Schicht- 
krystall) for certain cases of low misfits with growths of alkali halides 
upon each other. 

In the majority of cases which are studied by electron diffraction, 
growth occurs in the form of three-dimensional crystallites, independent 
of the misfit. This applies, for example, to the growth of alkali halides 
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upon silver, calcite and sodium nitrate (Schulz 1952 b,c). Metallic 
deposits from the vapour mostly grow in this way, although it has been 
found by Newman (results to be published) that lead and thallium deposited 
on to an atomically flat (111) silver surface at room temperature spread 
considerably (perhaps even as monolayers) during the growth of the first 
few angstréms in average thickness. This is illustrated by the diffraction 
pattern of fig. 17. The continuous streaks arise from the lead, or from 
secondary scattering involving lead and silver, and the few elongated 
spots from the silver substrate. This pattern should be compared with 
fig. 18, which shows the pattern obtained from the same thickness (2 A) 
of copper on a similar silver surface (Newman and Pashley 1955). Separate 
nuclei of copper have formed, giving rise to a pattern of distinct, although 
diffuse, elongated spots. It is estimated that these nuclei are in the form 
of plates which are about 124 high and about 60A across. Seifert (see 
Gomer and Smith 1953, p. 361) has expressed doubts concerning the 
electron diffraction evidence on nucleation effects; he suggests that a 
real monolayer exists in addition to the separated nuclei, and that this 
monolayer might remain undetected by electron diffraction. This is 
certainly not the case in most experiments, since any monolayer would 
give an observable streak pattern. 1n the case of copper on silver quoted 
above, it was found that separated nuclei were formed when the average 
thickness of the deposit was only about one-third of that of a single atomic 
layer, as determined by a radioactive tracer method. There was not 
enough copper present to form a monolayer of appreciable area, in 
addition to the nuclei. This is also true of the deposits of alkali halides 
studied by Schulz, although his thickness estimates were less accurate. 

The study of layers grown by methods other than from the vapour 
phase is difficult in the thickness range below 10 A, because of the lack 
of control which results from the inability to examine the layers during 
their growth. What evidence there is, suggests that separated nuclei 
are formed in most cases. Oriented chemically grown layers of silver 
bromide on silver form independent crystallites (see fig. 19), the size of 
which is estimated to be 404 high by 704 across in a particular case 
(Newman and Pashley 1955). There appear to be no examples of studies 
of very thin layers of metals electrodeposited on to atomically flat 
substrates. In most examples of oriented growth from solution, microscope 
observations reveal that large isolated crystallites form on the substrate, 
but the average thickness of the deposit is always high. It is not usually 
known whether any of the deposit in the form of monolayers or otherwise. 
remains between the crystallites which are visible under the microscope. 
For one or two examples of alkali halides grown from solution on to mica, 
specimens have been examined both by optical miscrocopy and by 
electron diffraction (Lisgarten 1954). The appearance of strong mica 
diffraction patterns shows that a large part of the mica surface has not 
been covered by the deposit, and that aggregation into large crystallites 
is complete. 
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5.2. The Mechanism of Formation of a Nucleus 


The occurrence of separated crystallites in a layer evaporated on to 
a flat substrate surface must result from the mobility of the deposit 
atoms. The manner in which this mobility leads to nucleation is not 
well understood, but from the various remarks which have been made 
in the literature, it may be said that there are two distinct possibilities. 
These may be termed as ‘gradual nucleation’ and ‘spontaneous 
nucleation ’. 

In a ‘gradual nucleation’ mechanism the deposit atoms are mobile 
over the substrate until they lose sufficient thermal energy and become 
trapped on the surface, either at a random lattice site or at some preferred 
place such as a hole, edge or step. Once trapped, the deposit atom acts 
as a preferential point for the trapping of other deposit atoms, so that a 
nucleus, which it is assumed can be oriented, is built up. By this 
mechanism both the number and size of the nuclei would increase 
with any increase in the deposit thickness. 

In a ‘spontaneous nucleation’ mechanism crystallization occurs via 
an intermediate stage, which could be one of an amorphous or liquid 
state. When the amount of deposit material is increased to a certain 
value a crystallization process occurs, and nuclei are spontaneously 
formed. Again it is assumed that the nucleus can be oriented. This 
type of process might form a large number of nuclei of almost the same 
size, so that for a certain period of the growth the number of the nuclei 
would increase, but their sizes would not increase appreciably. 

What evidence there is tends to favour some kind of ‘ spontaneous 
crystallization ’ mechanism in some cases. According to Schulz (1952 a), 
the nucleated deposits of alkali halides upon alkali halides first give 
observable diffraction patterns for average thicknesses as low as 2A. 
The linear size of the nuclei at this stage is already at least 100A. As 
the thickness of the layer is increased to about 100A, no marked change 
occurs in the size of the nuclei, but the increase in pattern intensity 
indicates a large increase in the number of the nuclei. Newman and 
Pashley (1955) find similar effects with copper deposited on to silver 
(111) surfaces. The diffraction pattern is observed continuously during 
the growth of the copper, which takes place at a rate of only about 
0-5A/min in average thickness. No visible change occurs on the 
diffraction pattern during the deposition of the first 0-84; a pattern of 
spots due to oriented copper nuclei of about 60A across then appears 
and these spots rapidly increase in intensity for further small increases 
in the deposit thickness. The number of nuclei increases, but their 
apparent size remains the same. This seems to be inconsistent with a 
‘gradual nucleation ’ mechanism, but could perhaps be reconciled with 
some kind of ‘spontaneous nucleation’ process. It is necessary to 
consider what form the deposit could take before appreciable nucleation 
occurs (i.e. at 0-8A) and why it gives rise to no observable diffraction 
pattern. Since very bright spot patterns are obtained from layers of 
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1-04 thickness, it would be expected that if the layers of 0-54 were of 
the same form and different only in amount, an observable pattern would 
be obtained from them. Only tentative suggestions can be offered at 
this stage. One possibility is that a large proportion of the first 0-84 
of the deposit becomes trapped in holes, cracks or grain boundaries, 
where it will be shielded from the electron beam. Alternatively, all 
the copper may be contributing to the diffraction pattern, but this 
contribution is either (a) spread over the background of the pattern, 
or (b) coincides with the silver pattern. (a) would result from an almost 
random distribution of the copper atoms; (b) would arise from regular 
arrays (perhaps two-dimensional) of the copper atoms, matching in form 
and spacing that of the silver substrate. In view of the results discussed 
in the next paragraph, (a) seems the more likely. 

There is further evidence that, for a given substrate temperature, 
nucleation will not occur until sufficient material is deposited to form 
nuclei larger than a certain tritical size. The solidification temperature 
of very small droplets is lower than the value for the bulk material. 
Turnbull (1950) has made measurements on the solidification temperature 
of droplets, 50 in diameter, of various metals and finds pronounced 
supercooling effects. Takagi (1954) has observed similar effects by 
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electron diffraction with thin films of lead, bismuth and tin on a substrate 
surface. The solidification temperature is found to decrease as the 
average thickness of the deposit decreases. The results, for the case 
of lead, are given in fig. 20. When the layer is cooled below the solidifi- 
cation temperature, the diffraction pattern changes from one due to an 
amorphous or liquid state to one due to a crystalline deposit, which is 
oriented in some cases. With a behaviour of this kind, a layer condensed 
on to a substrate at a given temperature would be in a liquid or amorphous 
state until a thickness is reached at which the solidification temperature 
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of the layer equals the temperature of the substrate. Nucleation would 
then occur. This kind of behaviour has, in fact, been observed by 
Newman (work to be published) for the growth of tin on a silver (111) 
face at room temperature, the nucleation occurring at an average thickness 
of about 104A, to give a partially oriented layer. 


5.3. Influence of the Substrate Surface Topography 


When an atomically flat substrate surface is used, the orientation of 
any deposit will be controlled, in general, by the type of lattice plane 
parallel to the surface. When the substrate surface is rough on an 
atomic scale, it will influence the orientation of any overgrowth in a 
more complicated manner. The simplest case occurs when the substrate 
is bounded by well-defined facets ; for the growth of silver halides on 
silver (Pashley 1952 a) and of oxides on cadmium and zine (Lucas 1952) 
it has been shown that the orientations which occur on such facets are 
the same as those which occur on flat surfaces of the same lattice type. 
In other cases, the interpretation of the observations is less obvious. 
It is possible that any form of surface irregularity, ranging from that 
which results from the termination of a dislocation line to larger features 
such as steps or cracks in the surface, might influence an overgrowth 
orientation. Such effects would be particularly marked if the growth 
were to occur preferentially at surface irregularities. Unfortunately, it 
is extremely difficult to determine experimentally whether any such 
surface features are of importance. Optical microscope evidence often 
shows that overgrowth crystallites are preferentially located at steps 
or other irregularities on the surface, as shown in fig. 21. However, 
there is usually no evidence to suggest that the irregularities have been 
responsible for the observed orientation. Thus, in fig. 21, a number of 
crystallites have formed along a curved step on the surface, but the 
orientation of the crystallites has not been influenced by the curvature 
of the step. Wilman (1951) has stated that rotational slip results show 
that the existence of angles or edges or cleavage steps is quite unnecessary 
for epitaxial orientation to occur. This deduction is far from convincing, 
since it is based upon data concerning the way in which two crystallites 
of the same substance are locked together, when they are misoriented 
with respect to one another by a small amount. This is, in the author’s 
opinion, irrelevant to any general rules of epitaxy; the latter involves 
the growth of one substance upon the surface of another. 

Although no direct evidence of the effect of small-scale irregularities 
is available, there are indications that such effects are operative in some 
cases. The oriented growth of silver bromide on variously prepared 
silver surfaces (Pashley 1952 a) provides some evidence. Although the 
main orientations observed were not attributed to surface irregularity 
effects, some subsidiary orientations were interpreted in this way. This 
applied, in particular, to an orientation in which a (211) plane of silver 
bromide is parallel to a (100) plane of silver. The evidence that this is 
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initiated by surface irregularities (probably steps) is the following : 
(1) the (211) orientation does not occur on atomically flat (100) silver 
surfaces, but does occur on electropolished (100) surfaces, which are 
known to contain irregularities ; (2) the (211) orientation occurs on (010) 
and (001) facets on an etched (100) surface, but not all equivalent 
orientations occur simultaneously; this can only be explained by 
asymmetries in the fine structure of the facet surfaces. It is interesting 
to note that this (211) orientation corresponds to a good fit only along 
the direction of the step at which it is supposed to form. 


5.4. Polymorphism in Thin Layers 


The measurement of the lattice spacings of small nuclei has been 
considered in § 4.3, it being shown that all evidence suggests that these 
spacings are nearly normal and are not constrained to fit that of the 
substrate. The structure of small oriented nuclei is not, however, always 
the same as that found for the bulk material. Rubidium bromide, 
which normally crystallizes with the rocksalt structure, is found to have 
the caesium chloride structure when deposited on to a silver substrate 
(Schulz 1951 a), the misfit being then zero; thallium chloride has a 
rocksalt structure instead of its normal caesium chloride structure when 
deposited on to potassium bromide (Pashley 1952 b). 

Schulz (1951 c) has studied deposits of caesium and thallium halides 
on various substrates. Although these salts normally oceur in bulk 
with the caesium chloride structure, the thin surface deposits often have 
the rocksalt structure. The complete results are summarized in table 20, 
which includes the values of the percentage misfits involved for both 
structures. The rocksalt type structures orient on the rocksalt type 
substrates in parallel orientation ; the caesium chloride type structures 
grow with their (110) planes parallel to the substrate (001) surface, 
and their [001] direction parallel to the [110] or [110] direction of the 
substrate. 

It is clearly shown that the substrate has some influence upon the type 
of structure which occurs in the deposit, but the effect does not seem 
to be associated with a small misfit. In some cases (e.g. CsI upon mica), 
the normal structure is observed although a smaller misfit (+-5°% instead 
of +25) would occur for the abnormal structure, which is observed 
upon other substrates. 

No other systematic evidence of the occurrence of abnormal structures 
in oriented layers has been reported. A close-packed hexagonal structure 
has been observed in an oriented copper layer electrodeposited on to a 
single crystal of B-brass (Takahashi 1952, 1953), the normal copper 
structure being face-centred cubic. This case corresponds to a zero 
misfit (see fig. 10), but it is possible that the transition is influenced by 
the presence of small amounts of some impurity (particularly zinc) in. 
the electrolytic bath. 
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5.5. Variations with Thickness 


An oriented layer will usually persist only up to a certain thickness, 
after which further deposition gives rise to less regular growth. In some 
cases the transition occurs when the average thickness is only about 104, 
while in other cases layers of several thousand angstroms are well oriented. 

Most of the evidence concerning thickness effects has been obtained 
for the growth of metal films, for which pronounced temperature effects 
are observed (see § 3.6). For metals (silver has been studied in most 
detail) deposited from the vapour phase on to various heated substrates 
(rocksalt, mica, etc.) three main effects can occur during thickness 
increase: (1) twinning; (2) changes in orientation; (3) changes in 
surface topography. 

Twinning on all four (111) planes occurs with most of the face-centred 
cubic metals. Where systematic tests have been carried out, it is found 
that the twinning does not occur in the thinnest deposits, but only 
develops as the thickness is increased above a certain value (Kirchner 
and Cramer 1938, Kehoe (results to be published)). One point of interest 
is that there is no conclusive evidence of the occurrence of second (or 
higher) twinning, i.e. twinning occurring on one set of (111) planes to 
give a new orientation (first-twin) which, at a later stage, twims on 
another set of (111) planes to produce a further orientation (second-twin). 
Such an effect produces sixteen orientations in addition to the initial or 
parent orientation. ‘ First-twinning’ gives rise to just four of these 
extra orientations. The orientations which arise may be readily deduced 
from the following expressions. A first-twin on the (111) plane has its 
(hkl) plane parallel to the (—h+2k+ 21, 2h—k-+-2l, 2h+2k—l) plane of 
the parent, and for a first-twin on a (111) plane the corresponding plane 
is the (—h+2k—2l, 2h—k—2l, —2h—2k—l). If a first-twin on the 
(111) plane is subsequently twinned on its (111) plane, the resulting 
second-twin has its (hk/) plane parallel to the (h—8k—4l, —8h-+-k—4l, 
. 4h-+-4k—'l) plane of the substrate ; there are corresponding expressions 
for other combinations of twin planes. 

Goswami (1954) claims to have observed second-twins with silver 
layers on a heated (001) rocksalt cleavage surface. When first-twinning . 
occurs, it is often accompanied by the appearance of orientations involving 
a (211) type silver plane parallel to the (001) rocksalt surface. Extra 
spots due to these orientations are visible on fig. 4 of a paper by Kirchner 
and Cramer (1938), and the same orientation is found by Kehoe (results 
to be published). Goswami points out that some of the second-twins 
have a (211) type plane nearly parallel to the surface and so concludes 
that these orientations are due to second-twinning. However, the 
application of the above formula shows that it is the (447) plane which 
is exactly parallel to the (001) surface. The (112) plane makes an angle 
of 33° with the surface. The observations are not consistent with a 
difference as large as this. It is therefore difficult to reconcile the (211) 
type orientation with second-twinning, although it does seem to have 
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some relation to the first-twinning, since it only appears when the twinning 
is present. In terms of the rocksalt substrate, this orientation is (211) Ag 
parallel to (001) NaCl, with the [017], [011], [111] or [111] Ag directions 
parallel to the NaCl [110]. The four positions are geometrically equivalent. 

There is evidence, in some cases, that the orientation of a deposit 
improves as it is built up. Randomly oriented material is formed during 
the early stages of growth, but the material deposited at a later stage 
is all well oriented. Rudiger (1937) has illustrated the effect quite 
clearly. Gold was deposited on to a thin sheet of mica at 380°c and 
examined by both reflection and transmission. The transmission pattern 
consisted mainly of rings, indicating that the bulk of the film had a fibrous 
structure, whereas the reflection pattern indicated complete orientation 
of the gold on its upper surface. Other evidence has not always been 
so decisive. When stripping methods are used as a means of study of 
thickness variations, spurious results can be obtained because the thin. 
layers (<200A) become randomly oriented as a result of the stripping 
process (Bruck 1936). Further, the dependence upon thickness can be 
influenced by rises in temperature which can occur during the deposition 
as a result of radiant heat from the source or as a result of the thermal 
energy carried by the molecules. Because of these factors, the behaviour 
of the deposits depends considerably upon the experimental conditions. 
If a slow rate of deposition is used, and the radiation effect is reduced to a 
minimum by the insertion of suitable shields, any temperature rises can. 
be made negligible. In some circumstances, the rises in surface tempera- 
ture can be as high as a few hundred degrees Centigrade (Wilman 1955). 

In view of the uncertainties involved, the majority of evidence con- 
cerning orientation improvements with thickness is not decisive. Even 
where the improvement is established as being a genuine effect of an 
increase in the thickness of the deposit, it is usually not possible to 
decide whether it has resulted from some kind of preferential growth 
on to the better oriented nuclei, or whether a complete recrystallization 
of the film has occurred at a certain stage of the deposition. More 
experiments on this aspect of oriented growth are necessary before any 
further conclusions can be drawn. 

Since, in most cases, the metal is deposited in the form of separated 
crystallites, the initial layers are rough on an atomic scale. As the 
crystallites grow together to form a coherent film, there is often a 
tendency to form a fairly smooth surface. Silver on heated rocksalt 
behaves in this manner; in extreme cases, an atomically flat surface 
may be formed (Tull 1951). The author has studied in detail (results 
to be published) the growth of silver on mica heated to about 250°c. 
With a rate of deposition of a few hundred 4/min the deposit starts off 
as partly oriented in (111) double positioning ([110] and [110] Ag parallel 
to mica [010]), and partly random. This is followed by twinning on all 
(111) planes and the disappearance of random orientation in the surface 
layers. As the thickness reaches several hundred angstroms, the twinning 


228 D. W. Pashley on the 


becomes less pronounced and the surface starts to become smooth. The 
surface orientation gradually changes from double to single positioning 
({110] Ag only parallel to mica [010]). With a thickness of more than 
about 10004 the surface is atomically smooth and it is oriented with 


complete single positioning. } 


5.6. Perfection of Orientation 


The perfection of orientation varies considerably from case to case. 
Sometimes the deposits give rise to diffraction spots which are appreciably 
arced, indicating that the individual crystallites are aligned to an accuracy 
of only a few degrees. In other cases, the bright sharp Kikuchi line 
patterns from the deposit show that the alignment is much more accurate. 
Figure 22 shows a high angle Kikuchi pattern from a layer of silver on 
mica. The sharpness of the lines indicates that the variation in alignment 
of the silver crystallites is less than about one-fifth of a degree. 

Some interesting misorientation effects are sometimes observed. It is 
found that a layer studied by the reflection technique gives a pattern of 
fairly sharp spots at one azimuth, but appreciably arced spots at the 
perpendicular azimuth (cf. figs. 23 and 24). Clearly, a misorientation 
has occurred as a result of a partial rotatory degree of freedom about 
one preferred axis parallel to the surface. This effect has been observed 
by Evans and Wilman (1950), with a layer of zine oxide obtained by 
heating a cleavage (110) surface of zinc blende. The oxide is oriented 
with its (1013) plane parallel to the surface. When the electron beam is 
sent along the oxide [010] direction, the diffraction spots are prominently 
arced ; along the perpendicular azimuth, the diffraction spots are com- 
paratively sharp. Their interpretation is that a disorientation of initially 
well aligned nuclei occurs as a result of rotational slip on the (1210) 
planes which are perpendicular to the [010] direction, so that an 
‘epitaxial strain’ is relieved. This strain is supposed to result from 
the initial layers of the deposit having a pseudomorphic structure. 
There is no other evidence to support this hypothesis of the pseudo- 
morphism of the zine oxide layer. The explanation of Evans and Wilman 
does not seem very likely, since a similar effect has been observed 
(Pashley 1952 a) in a case where the misfit is zero and no epitaxial strain 
can possibly exist. This example is illustrated in figs. 23 and 24, which 
are patterns obtained from an oriented silver bromide layer on a (110) 
surface of silver. Along the [110] AgBr direction the diffraction spots 
are prominently arced (fig. 23); at the [001] azimuth, the spots are 
sharp. Misorientation occurs about the [110] axis. 

The origin and significance of these misorientation effects are still, 
therefore, a matter of speculation. 


§ 6. THEORIES OF EPITAXY 
All theoretical considerations of the occurrence of epitaxy have been 


based upon the concept of the geometrical fitting between the lattices 
of the substrate and of the overgrowth. As a result of his extensive 


Study of Epitaxy in Thin Surface Films 229 


studies, Royer (1928) put forward the three following rules for epitaxy : 
(1) That there exist lattice planes with networks, elementary or multiple, 
which are identical in form and of nearly the same dimension in the two 
structures. It is implied that these two networks grow parallel to one 
another. (2) That, where ionic crystals are involved, the ions of the 
overgrowth should take up positions which corresponding ions of the 
substrate of the same polarity would have occupied had the substrate 
continued to grow. (3) That the substrate and overgrowth crystals 
should have the same type of bonding. 

_ The extensive data on epitaxy which has been accumulated since 
Royer’s early work have shown quite definitely that the third rule is 
unnecessary. Thus, for example, ionic and metallic crystals can be 
mutually oriented (see § 3.4). In accord with the later views of Royer 
(1954), this rule must be rejected, at least in its original form. 

The oriented growth of various metals on heated rocksalt (see § 3.6) 
could not be explained directly in terms of Royer’s rules, since the 
orientation predicted by rule (1) did not occur. Bruck (1936) attempted 
to explain the results by saying that the orientation is such as to make 
the sum of the distances of the substrate ions from the deposit atoms a 
minimum when considered over one unit cell only. This rule has been 
criticized by Thomson and Cochrane (1939), who point out that what 
would be a favourable orientation for the first cell might be an unfavour- 
able one for nuclei of greater lateral extent. Also, the effect of both 
positive and negative ions being present was not considered, and Engel 
(1953) has attempted to develop the rule to take this into account. 

Dankovy (1946, see also Pinsker 1953, p. 220) has attempted to establish 
the condition for the growth of an oriented crystal nucleus, strained to 
fit the substrate, by considering the energy of deformation which results 
from the misfit. However, this calculation cannot be carried out in 
most cases, because of the absence of the necessary lattice data. 
A calculation based upon the hypothetical growth of sodium chloride 
on sodium, under the action of chlorine gas, was carried out, but this 
led to no satisfactory application to other cases. 


6.1. The Menzer Mechanism 


Some attempts at explanation of the occurrence of large misfits have 
been based upon the idea that an oriented layer corresponding to a good 
fit occurs during the initial stages of growth and that subsequent growth 
gives rise to different orientations. The most notable of these mechanisms 
is that put forward by Menzer (1938 a, b, ¢) for the growth of nickel 
and silver on rocksalt. This is based upon the data of Bruck (1936), 
which showed that such layers are prominently twinned on their (111) 
planes when the film thickness is several hundred angstroms. Since 
the twinned lattices have (221) planes parallel to the rocksalt surface, 
Menzer suggested that the initial nuclei of the deposit grow with (221) 
planes parallel to the surface of the rocksalt, there being four geometrically 
equivalent positions for these nuclei (fig. 25). As deposition continues, 
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these nuclei are supposed to twin on their (111) planes and grow together 
to form an interpenetrating lattice which eventually develops into the 
parallel orientation. Goche and Wilman (1939) considered this mechanism 
to be consistent with their observations on silver films on rocksalt. The 
reciprocal lattice corresponding to the observed pattern consists of a 
body-centred cubic array of points (see fig. 6 of Goche and Wilman’s 
paper) due to the parallel orientation, plus extra spots along the body 
diagonals ; these spots divide the diagonals into six equal parts, Only 
some of the extra spots arise directly from the (221) or twin orientations ; 
Menzer explained the remainder as due to the junction region between 
the parallel and twin lattices in which there is an incomplete filling of 
the (111) planes. This region is thus supposed to have a periodicity 
which is three times that of the normal lattice. The intensity of the 
extra spots indicated that the junction regions formed an appreciable 


part of the film. 
Fig. 25 


NoC2 [100] 


Nace [o1o} 

A plan view of four (221) silver nuclei on an (001) rocksalt surface. The 
boundary faces, all of (111) type, make the following angles with the 
rocksalt surface: (111), 15°54’; (111), 54°45’; (171) and (111), 
78° 64’. 


According to this mechanism, the misfit of the initial nuclei on the 
rocksalt substrate is only —7°% for nickel and + 9° for silver (see fig. 26). 
However, it should be pointed out that these misfits apply only if one 
silver (or nickel) distance is compared with two sodium distances along 
the silver (or nickel) [114] and if three silver (or nickel) distances are 
compared with two sodium distances along the silver (or nickel) [110]. 
This is really no better a fit than occurs for the parallel orientation, 
which is given in fig. 27. Along all of the principal directions the misfit 
can be regarded as +-9°,, for silver (or —7°% for nickel) if three silver 
(or nickel) distances are compared with two sodium distances. 


Study of Epitaxy in Thin Surface Films EH | 


It should also be pointed out that the mechanism only accounts for 
the twinning of the initial nuclei on their (111) planes. Twinning on the 
other (111) type planes, particularly the (111) (see fig. 25), of the (221) 
nuclei would give rise to orientations other than the parallel one: these 
are not observed. : 


7 >) D oa 
Fig. 26 Fig. 27 
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SURFACE , 
The fit of the (221) silver plane The fit of the (001) silver plane on 
(Menzer hypothesis) on a a rocksalt (001) surface. 


rocksalt (001) surface. 


6.2. The Engel Theory 


Engel (1952, 1953) has put forward a theory for the orientation of 
metals upon salt substrates, in which she attempts to explain both the 
observed fits and the dependence of the orientation upon the substrate 
temperature. This dependence (Bruck 1936) amounts to there being 
a minimum substrate temperature (the epitaxial temperature) for the 
occurrence of complete orientation. It is suggested that the first stage 
in the formation of an oriented film is that the atoms in the first layer 
become ionized, so that a two-dimensional chemical reaction occurs. 
A monolayer of the deposit—metal salt then exists at the interface between 
the substrate and the overgrowth. According to this theory, the 
occurrence of orientation depends upon the possibility of suitable 
ionization processes and upon the misfit of the intermediate salt with 
both the substrate and the deposit. An empirical approach indicates 
that an approximately linear relation exists between Bruck’s observed 
epitaxial temperatures for various metals, and the ionization potential 
corresponding to the most commonly observed ionization state of the 
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metal. Engel suggests that this can be used to predict the epitaxial 
temperature for other metals. ; 

As an illustration of the fits involved in this mechanism, the case of 
silver on rocksalt is quoted. The misfit of the silver chloride on rocksalt 
is only —2°, and the misfit of the silver on the silver chloride is again 
only —2°%, if three cell lengths of the silver chloride lattice are compared 
to four cell lengths of the silver. However, this does not amount to any 
better fit than occurs without the silver chloride layer, since the 
application of this latter fitting procedure gives a misfit of only —3%, 
for silver direct on to rocksalt. 


6.3. The Frank and van der Merwe Theory 


Van der Merwe (1949) and Frank and van der Merwe. (1949 a, b, c) 
have put forward a theory of orientation, in which it is assumed that 
the initial stage of growth of an oriented deposit is the formation of an 
immobile monolayer of regular atomic pattern. This monolayer is 
considered to be homogeneously deformed to fit on to the substrate 
(i.e. have the same spacing), in order to form a suitable embryo. The 
theory is developed so as to lead to predictions regarding the necessary 
conditions under which such an embryo can be formed. A one-dimensional 
dislocation model is used, but reasons for justifying the extension of the 
results to two dimensions are given. The overgrowth is represented by 
, a row of identical balls connected by a row of identical springs (force 
constant .); the substrate is represented by a periodic force acting on 
the balls, the first harmonic term (4W) of which is taken to represent 
the corresponding potential energy. It is shown that the lowest energy 
state of the system is the one with a monolayer constrained to fit the 
substrate, provided the misfit is less than a certain critical value given 
by (6/4—1)ctitiea=2/7lo¢ Where a and b are the natural spacings of the 
substrate and deposit respectively and J,2=a?/2W. This critical misfit 
should be about 9% in an average case, but would vary widely according 
to the relative values of «1 and W in a particular case. Above this critical 
misfit, dislocations will occur in the monolayer, but there remains an 
activation energy for their generation, which will be above zero for 
misfits less than a second critical misfit value of 1//,, which will be about 
14°, in an average case. Hence, below this second critical misfit, it is 
possible to form a metastable monolayer at sufficiently low temperatures. 

Once this oriented monolayer has formed, a macroscopically thick 
oriented film can grow by repetition of the process. It was considered 
that this would give rise to an initially pseudomorphic film, which would 
undergo some transition process at a certain stage to give a strain free 
bulk deposit which remains oriented. 

To summarize, according to this theory a necessary criterion for 
orientation is the formation of an initial pseudomorphic monolayer 
over the substrate surface. This will occur provided the misfit is less 
than a certain critical value, which is about 14% in an average case. 
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6.4. Discussion 


Since all of the theories involve the concept of good fitting between 
the substrate and overgrowth, the evidence concerning this will now be 
discussed. 

Although many authors have explained the occurrence of a particular 
orientation as due to the existence of a small misfit, an examination of 
tables 1 to 18 shows that many cases of epitaxy involving large misfits 
are found. Any statistical analysis of all of the results is difficult, if not 
pointless. One of the best ways of considering the question is to examine 
cases where a series of similar deposits, with different spacings, has been 
studied on the same substrate. The results when both the substrates and 
the deposits are alkali halides are given in table 1. It is seen that although 
Royer and Sloat and Menzies found fairly well-defined misfit limits 
(e.g. —138% to +9% for a KCl substrate) for growth from aqueous 
solution, Schulz and Ludemann found no such limits for growth from 
the vapour phase. In these latter experiments, any alkali halide is 
found to orient on any other alkali halide and the complete results of 
these two authors indicate that epitaxy occurs for misfits covering the 
range —43°% to +90°%. These are the limits set by the lattice constants 
available. Similar results followed from observations of the growth of 
alkali halides upon mica. Furthermore, it has been pointed out on 
several occasions (e.g. Thomson 1948) that the observed orientation 
is not always that which corresponds to the best fit. This applies 
particularly to the growth of metals upon alkali halides. The orientations 
corresponding to good fits are not observed, although they are observed 
when alkali halides are grown upon the metals (Johnson 1950 b, c, d, 
Schulz 1952 c). 


-80 -60 -40 -20 ° 420 +40 +60 +80 +100 +120 
NON-MULTIPLE SMISFIT 


The regions of good fit introduced by the consideration of multiple fits. 


In view of these conclusions, the concept of the need for a small misfit 
must be qualified considerably. The most obvious and most important 
conclusion is that the occurrence of a small misfit is not an essential 
condition for the appearance of orientation. Thus Royer’s first law is 
invalid,.at least in its original form. 

When an orientation with a large misfit is observed, it is common 
practice to consider multiple fits, ie. the fit of m atom distances in the 
deposit with n atom distances in the parallel direction in the substrate. 
This is, however, of dubious value, because a certain flexibility of 
application soon leads to an absurdity. If one allows only integral 
values of m and n of three or less, a misfit tolerance of +12% would 
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allow almost all non-multiple misfits between —70% and +130% to 
be classed as good fits, as illustrated in fig. 28. The only appreciable gap 
is in the region of 20°, misfits and this would immediately be filled if 
m and n values of four were allowed. The significance of these remarks 
is illustrated by the fact that values as high as seven (Royer 1932) are 
sometimes put forward. Thus arguments concerning multiple fits must 
be treated with caution. 

The Menzer mechanism (see § 6.1), which represents another attempt 
at explaining the occurrence of large misfits, is based upon the experiments 
of Lassen and Bruck (1935 a). In their work, silver films of several 
hundred angstréms in thickness were stripped from rocksalt substrates 
and examined by transmission ; parallel orientation plus twinning was 
found. There was no direct evidence of the mode of growth of the 
initial nuclei. The growth of very thin layers of silver was studied by 
Kirchner and Cramer (1938), who used the reflection technique. They 
showed that twinning is not very prominent during the initial stags, 
but that it only becomes pronounced after the growth of the first 204 
in average thickness. Similar results were found by Uyeda (1942). 
Later, Raether (1946, 1951) and Tull (1951) stripped layers of silver 
from rocksalt and examined, by the reflection technique, the silver 
which had been in contact with the rocksalt. Only silver in parallel 
orientation was observed. Raether (1951) concluded that Menzer’s 
mechanism does not apply. Recently in this laboratory, Kehoe (results 
to be published) has examined the growth process in detail, by con- 
tinuously studying the formation of deposits in situ in the diffraction 
camera. He finds that in the thinnest observable layers the silver is 
in parallel orientation and that twinning occurs only at a later stage. 
No (221) oriented nuclei occur as an initial stage of growth. These 
results are illustrated by figs. 29 to 31, which are from the same specimen 
at different stages of growth. In fig. 29, which corresponds to an average 
thickness of silver of only about 24, the elongated spots are from the 
rocksalt substrate and the diffuse spots forming a centred rectangular 
pattern are from the silver in parallel orientation. Twin spots begin 
to appear at a thickness of about 104 and are clearly visible on fig. 30, 
which is from a layer of about 204 in thickness. This pattern is due 
entirely to the silver. At much greater thicknesses the silver surface 
becomes smoother, as indicated by the elongation of the main spots on 
fig. 31. 

Further evidence concerning the Menzer hypothesis was obtained by 
Kehoe. He has studied growth of silver, copper and gold on several 
substrates and finds parallel orientation in all cases, at sufficiently high 
temperatures. Some of these cases (see table 21) would involve quite 
large misfits if they were interpreted in terms of the Menzer mechanism. 

These results show, beyond all doubt, that the Menzer mechanism 
does not apply to the growth of a number of metals (including silver) 
upon alkali halides. The initial stage of growth is the formation of 
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nuclei in parallel orientation and twinning only occurs at a later stage 
of growth. 

As a result of the data presented in tables 1 to 18, and on the basis of 
the above arguments, it is concluded that : 


(1) A small misfit is not an essential condition for the occurrence of an 
oriented overgrowth. 


(2) The orientation which occurs is not always that which corresponds 
to the best geometrical fit. 

(3) Many cases of observed large misfits are not satisfactorily explained 
by the assumption of an intermediate orientation of smaller misfit, 
and the occurrence of very large misfits must be taken into account 
in any general theoretical treatment of epitaxy. 


It follows from these conclusions that the Frank and van der Merwe 
theory (see § 6.3) cannot hold in general, since it predicts certain critical 
misfits for the occurrence of orientation. It is, however, necessary to 
consider whether this theory has a more limited application. Smollett 
and Blackman (1951) have criticized the theory on the basis that it 
predicts constrained monolayers which would not necessarily be stable. 
Certain examples of alkali halide monolayers in parallel orientation on 


Table 21. Misfits for Observed Cases of Metals Oriented 
upon Alkali Halides (Kehoe, results to be published) 


°, misfit % misfit 
Substrate Deposit for the according to the 
parallel orientation Menzer mechanism 
NaCl Cu —36 —4 
NaCl Au —28 +8 
NaCl Ag —28 +9 
KBr Cu —45 —18 
KBr Ag —38 —7 
KCl bit —42 —14 
KI Cu —d0 =23 


an alkali halide (100) surface were considered in detail. The effective 
elastic constants of the monolayers were found; it was deduced that 
in some cases where oriented growth is known to occur, a monolayer 
constrained to fit on to the substrate would not be dynamically stable. 
It was therefore concluded that the orientation observed experimentally 
does not arise from an initial uniformly strained monolayer (for further 
comments on this point see Frank 1951 and Blackman 1951). Further, 
again with reference to the growth of alkali halides upon each other 
(see table 1), it is found that critical misfits are observed only in the 
case of growth from solution. However, it would be expected that the 
Frank and van der Merwe theory would apply to the vapour growth, 
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rather than growth from solution, since the theory did not take into 
account any possible influence of the solvent upon the growth. The 
results of Sloat and Menzies (1931) suggest that the presence of a solvent 
tends to reduce the chance of orientation rather than to increase it. 
Thus, in the case where the theory would be expected to be most 
applicable, there is no agreement with experiment. As is suggested by 
the calculations of Smollett and Blackman (1951), the discrepancy is 
probably due to the fact that the first stage of orientation is not the 
formation of a strained monolayer. The latest experimental evidence 
(see § 5.1) shows that in many cases, including some examples of growth 
of alkali halides upon each other, monolayers are not formed at all, but 
that surface migration of the deposit atoms leads to growth of small 
oriented three-dimensional nuclei. The Frank and van der Merwe 
theory would then not apply because it is based upon an initial mechanism 
which is not relevant to these cases. 

The experiments of Kehoe which are quoted above also throw doubt 
upon Engel’s theory (see § 6.2) of the effect of substrate temperature 
on the orientation of metals upon salts. Kehoe finds that good orientation 
occurs at temperatures appreciably lower than Bruck’s (1936) epitaxial 
temperatures. \Precautions were taken to ensure that no appreciable 
heating occurred during the deposition of the layers. The rates of 
deposition were lower than those used by Bruck, and smaller thicknesses 
were examined. The results suggest that the epitaxial temperatures 
observed by Bruck were a characteristic of the ability of the metal to- 
orient upon itself, after the formation of initial oriented nuclei, rather 
than to orient upon the substrate. This is presumably controlled by an 
effect associated with surface mobility. 


§7. SUMMARY AND Discussion 


The electron diffraction method has contributed in two main ways to 
the subject of epitaxy. It has supplied many new examples of oriented 
growth and has provided certain information about the mode of growth 
of the initial nuclei of an oriented deposit. Although it has been 
concluded that a small misfit is not an essential criterion for the 
occurrence of epitaxy, it is not intended to imply that the misfit value 
has no significance. The results of Royer and of Sloat and Menzies 
(see tables 1, 2 and 3) provide very convincing evidence that the misfit 
value is significant under certain conditions. However, the writer is of 
the opinion that the concept of small misfit values has dominated the 
subject more than is justified by the experimental facts. Other approaches 
to the problem are needed. 

The concept of basal plane pseudomorphism, which is linked with the 
small misfit criterion, has likewise been applied too widely, in view of 
the limited experimental backing. As is argued in § 4, none of the 
original evidence is valid, and more recent experiments have failed to 
supply any support for the concept. 
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The failure of existing theories in supplying a satisfactory explanation 
of epitaxy is that they are based upon unsuitable models for the growth 
of the initial nuclei of the deposit. In the case of the Menzer mechanism, 
the explanation was based upon effects (e.g. twinning) which are now 
known to occur only after the initial oriented nuclei have been formed. 
In the case of the Frank and van der Merwe theory, the concept of 
pseudomorphic monolayers was used. These monolayers certainly do 
not form in many cases. It therefore follows that the most important 
aspect of the study of epitaxy is the consideration of the mode of 
formation of the initial oriented nuclei. This is likely to vary con- 
siderably, according to the method used for the preparation of the 
deposit, as well as to other experimental conditions. Chemically grown 
deposits require special attention, because the substrate undergoes 
changes during the growth of a surface film. 

Both experimental and theoretical studies of the nucleation problem 
are highly desirable. As a first step, it is necessary to determine the 
size and form of the smallest stable nuclei, under a given set of growth 
conditions. The energy of all possible arrangements of these nuclei 
must then be considered, so that the most probable orientation. can be 
determined. This is a formidable task, and very little progress can be 
expected in the near future, because of the considerable lack of knowledge 
about the surface forces which are involved. However, from this point 
of view it might reasonably be expected that some form of preferred 
orientation will always occur provided the activation energies involved 
are not too high. Accordingly, as has already been pointed out by 
Schulz (see Gomer and Smith, 1953 p. 374), the actual problem should 
be to account for the failures to observe orientation. 

Many failures may be due to the limitations of the method of 
observation. In particular, the study of extremely thin layers has shown 
that good orientation sometimes occurs only during the growth of the 
first hundred angstroms (e.g. Schulz 1952a, for alkali halides upon 
themselves), or even just the first ten angstroms (e.g. Newman and 
Pashley 1955, for copper on silver). It is interesting to note that in 
both of these quoted examples, the poor orientation of the deposit 
occurs at a stage when the initial nuclei have just about covered the 
substrate surface. It is at this stage that the growth ceases to take 
place upon the flat substrate surface and starts to occur predominently 
on the existing deposit nuclei. Failure to observe orientation in these 
thicker films is therefore not directly related to growth of the deposit 
upon the substrate. ; , 
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§ 1. INTRODUCTION 


RECENTLY there has been shown an increased interest in the thermo- 
electric properties of semiconductors, partly because of possible technical 
applications of such new semiconducting materials as InSb (Ross and 
Saker 1955, Saker, Cunnell and Edmond 1955) and Bi,Te, (Goldsmid and 
Douglas 1954, Goldsmid 1955). It is well known that the theory of these 
effects is very complicated, even if one considers only the electronic 
component (compare a forthcoming review article by V. A. Johnson). 
However, it is impossible to account for the observed steep increase of 
thermoelectric power at low temperatures, if one considers only this 
component. In the present paper we are therefore considering the lattice 
component of the thermoelectric power and at the same time the lattice 
component of the thermal conductivity which follows from the same 
equations. In §2 we discuss the various approximations which are 
usually made in a discussion of transport properties and from that dis- 
cussion it is clear that it is well nigh impossible not to make some 
assumptions. We therefore decided to investigate as comprehensively as 
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possible the various scattering mechanisms, but to neglect the influence 
of deviations from spherical energy surfaces and many-electron phenomena. 
We also restrict our discussion to the case where only one type of carrier 
is present, although it would be straightforward to extend the theory to 
the case of both hole and electron conduction. 

One must consider simultaneously phonon scattering and electron 
scattering. Following Makinson (1938) we assume electron—phonon 
scattering to be active at high and electron-impurity scattering to be 
active at low temperatures, while at low temperatures we consider 
scattering of phonons by boundaries, electrons, and impurities and at 
high temperatures by electrons, impurities and phonons. This gives us 
six temperature ranges to consider for which we give expressions for the 
thermal conductivity and thermoelectric power (see equations (168) to 
(179)). Most of the paper is concerned with the derivation of these 
equations, but the last section is devoted to a discussion of the results 
obtained and a comparison of these with the scarce experimental material 
available. It is shown that although qualitative agreement exists between 
theory and experiment, this is not true quantitatively. 


§ 2. GENERAL CONSIDERATIONS. 


* The quantum theory of solids has the reputation of being rather less 
respectable than other branches of modern theoretical physics. There is 
some justification for this view because the dynamics of many-body 
systems, with which it is concerned, cannot be handled without the use of 
simplifications, or without approximations which often omit essential 
features of the problem.” These opening sentences of Peierls’ recent 
monograph (1955) indicate the problems which confront us when we try to 
explain theoretically the properties of solids. Although the Drude- 
Lorentz theory (Lorentz 1909) as modified by Sommerfeld (1928) and 
Sommerfeld and Bethe (1933) can account in many cases qualitatively for 
the observed transport properties of metals and semiconductors, this so- 
called elementary electron theory (for an account we refer to Wilson 1953 
(in the following referred to as TM) Ch. I, Seitz 1940, Ch. IV, ter Haar 
1954, Ch. X) cannot quantitatively account for the great diversity of 
observed transport phenomena. 

For our further, discussion it is useful to list various assumptions intro-— 
duced in this theory. They are the following ones : 


(i) It is possible to define a relaxation time, 7, or a mean free path, /, 
for the charge carriers ; 

(ii) one may neglect the influence of the lattice, apart from the fact 
that it provides a mechanism by which a finite relaxation time may be 
realized ; 

(ili) if fy denotes the equilibrium carrier distribution function and f the 
carrier distribution function in the presence of external fields and/or 
temperature gradients, |(f—f )/fo|<1 ; 
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(iv) one may treat the carriers as if they were independent, that is, one 
may use the so-called one-electron approximation ; 


(v) the equilibrium distribution function f, depends on the carrier wave 
vector k only through the energy E ; 


__ (vi) the energy surfaces are spherical, that, # depends on k only through 
its absolute magnitude k. 


In the development of the elementary electron theory further assump- 
tions are introduced, such ‘as the degeneracy of the gas of conduction — 
electrons in the case of metals or the ‘ quasi-free ’ approximation (see 
eqn. (17) below). It is interesting to note that assumptions (i) and (iii) 
together suffice to prove (compare TM p. 196, Peierls 1955, § 6.1) that one 
may use for f the usual Lorentz expression (1905, 1909 ; see eqn. (15) 
below). 

Several attempts have been made to consider the influence of removing 
one or more of the restrictive assumptions, but to remove all, and to take 
into account as well the recently acquired knowledge about the macro- 
scopic structure of solids, is clearly a Herculean task and it seems that for 
some time to come one must be satisfied to tackle the problem piecemeal 
and to investigate the influence of the various factors separately, hoping 
that cross effects will not seriously upset the conclusions reached by this 
approach and that one’s approximations are at least consistent. That 
one can easily be led astray was illustrated by Sondheimer (1956) who 
showed that the work by the present authors (ter Haar and Neaves 1955) 
on the thermoelectric power of metals suffered from inconsistency in the 
approximations used. 

Various Russian authors have in recent years tried to eliminate 
assumption (iv), but up to now their theories do not seem to be in a state 
where quantitative numerical predictions can be made. We refer to survey 
papers by Anselm (1951) and Wonssowski (1954) for a discussion of these 
theories, and only mention here that it looks from their work as if this 
particular assumption does not affect the final results to any great extent 
(see, however, Bontsch—Brujewitsch 1955 and Pines 1955). One should 
have expected this, since in the case of metals the Pauli principle will 
effectively restrict the possible interactions between the conduction 
electrons, while in semiconductors the carrier density will be small 
enough for two-carrier processes to be negligible. 

An interesting point is in how far assumption (vi), especially when 
taken in conjunction with the ‘ quasi-free’’ approximation, .is justified. 
It is perhaps surprising that nobody seems to have attempted to introduce 
the vast amount of results from band theory into a systematic discussion 
of transport phenomena, apart from introducing an effective mass into 
the discussion (compare Peckar 1946, James 1949, Slater 1949) and from 
using a two-band approximation (e.g., TM p. 198, Hung 1950). Recent 
work on cyclotron resonance in germanium (Dexter et al. 1954, Lax et al. 
1954) has shown the necessity to introduce warped energy surfaces and the 
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so-called many-valley model is now being investigated (Zeiger 1959, 
Herring 1955, Zeiger, Lax and Dexter, private communication, Lax and 
Mavroides, 1955). 

The main interest of many recent investigations has been in the dis- 
cussion of assumptions (i) and (ii) and in the calculation of the relaxation 
time +. We refer to TM for an account of this work. The present paper 
falls into this category and we shall throughout introduce assumptions (iii) 
to (vi) as well as the ‘ quasi-free ’ approximation. Our investigations were 
instigated by the experimental results of Frederikse (1953) on the thermo- 
electric power of germanium at low temperatures. His work was an 
extension of the measurements by Lark-Horovitz et al. (1946) and by 
Middleton and Scanlon (1953). The latter experiments could readily 
be explained by using the elementary electron theory (Johnson and 
Lark-Horovitz 1953, see also Lautz 1953), taking into account that there 
is both electron and hole conduction and that degeneracy (Johnson and 
Lark-Horovitz 1947) may occur. However, Frederikse’s results could not 
be explained in that way and he suggested that there might be a contri- 
bution from the lattice current, the so-called Gurevich effect (1945); the 
same suggestion has been made independently by Herring (1953, 1954 b ; 
see also Klemens 1954, Parrott 1954) to explain similar experimental 
results obtained by Geballe and Hull (1954) and Geballe (1953). In the 
present paper we discuss the influence of lattice heat current and lattice 
non-equilibrium on both the thermoelectric power and the thermal 
conductivity. As in an earlier paper (ter Haar and Neaves 1955, to be 
referred to as A) we shall use Makinson’s method (1938, to be referred to as 
M) which was an extension of a paper by Wilson (1937) on the second order 
electrical effects in metals. It will be seen that, indeed, the influence of 
the lattice is of great importance, especially at low temperatures. A 
preliminary account of the results of some of our calculations on the 
thermal conductivity was given at the 1954 Amsterdam Conference 
(ter Haar and Neaves 1954). 

We should at this point make a few remarks about our earlier paper on 
the thermoelectric power of metals at temperatures around the Debye 
temperature. In that paper we calculated the Thomson coefficient y 
which is related to the thermoelectric power Q by the Kelvin relation 


Q=| WiTVaT Gen ea 


In our calculations we assumed (A, eqn. (25)) that we could neglect the 
electron non-equilibrium in calculating the lattice distribution function. 
However, Sondheimer (1956) has pointed out that this assumption leads to 
a violation of the Onsager relations and that one can thus no longer expect 
that the expression for Q derived from eqn. (1) will be the same as the one 
derived directly from the current density J in the presence of an electric 
field & and a temperature gradient d7'/dx, 


eJ =o|(e& +-dl/dx)— eQ dT/da], emia, haiti og VU)" 
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where ¢ is the Fermi-level, o the electrical conductivity and e the charge of 
the carriers. As a result most of the results of A will be incorrect. More- 
over, we should like to express our regret for not having discussed in A, 
Gurevich’s paper, the results of which in many respects agree with ours 

at least qualitatively. 

_ For an understanding of the physical background of the thermoelectric 
effect it is illuminating to consider briefly the phenomenological theory 
of these effects (TM §8.4, ter Haar 1954, App. II, de Groot 1951, 
Samoilowitsch and Korenblit 1954). Consider a piece of wire under the 
mfluence of an electric field € and a temperature gradient d7'/dxz. There 
will be both an electrical current density J and a heat current density w 
and we have the relations 


J=A(e€+dt/dx)+BdT/dx, ..... (3) 
w=O(e€+dt/dx)+DdTidx, .... . (A) 


where the gradient of the Fermi-level acts as an electrical field, since the 
Fermi-level is, indeed, the electrochemical potential (compare the dis- 
cussion on p. 229 of ter Haar 1954). 
If the experimental set-up is such that J=0 we get the thermoelectric 
power @ from eqn. (3), 
Q=—B/Ae. emer ge 64 a (5) 
If, on the other hand, d7'/dx=0, we are dealing with the Peltier effect, 
and the Peltier heat J7, which is the flow of heat per unit electric current 
is given by the equation 
T==ClA. ce ee en ee. (6) 


Finally, the Thomson coefficient y is defined as the negative of the 
coefficient of J(d7/dx) in the expression for the rate at which heat 
accumulates, dH/dt, in a unit volume of the wire. Now, we have 


dH/dt=J & —du/dx, eae 
or eliminating & from eqns. (3), (4) and (7) 
dH /dt=J?/A e—(J/e)(dl/dx)—(B/A e)J (dT'/dx)—J[d(C/A)/dT (dT /dz) 
+ (d/dx)[{(BC—AD)/A }(dT/dz)], etcs ee eS) 
or ya BiAe+d(C/A)jdT. 7. . ~. 2 «+ (9) 


From eqns. (5), (6) and (9) we see that, if the Onsager relations. are 


satisfied, that is, if 
BT=—Ce, weeny te ( 10) 


eqn. (1) and the other Kelvin relation 
TL ee ee ee ee 


are realized, and one can thus use any of the three eqns. (5), (6) or (9) to 
calculate the thermoelectric coefficients. 
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Before we discuss the results obtained by using Makinson’s equations we 
want to remind the reader that one can obtain from simple kinetic argu- 
ments (ter Haar 1956) the following expressions for the thermal conduc- 
tivity « and for the thermoelectric power Q, 


Ke=cvl’, ay taney ee need eee 


eQ=—f(e/n)(U)))..) a+ le ee 
where in both equations numerical factors have been omitted, Also ¢ is 
the specific heat per unit volume which is the electronic specific heat when 
we are dealing with pure electronic effects and which is the lattice specific 
heat when the lattice heat current dominates, /’ is the mean free path 
of the heat carriers (electron mean free path or phonon mean free path) and 
v their mean velocity. In eqn. (23) m is the number density of the charge 
carriers, / the electron mean free path, and f a factor indicating the 
efficiency with which the phonons carry the electrons along when we are 
dealing the contribution to Q due to the lattice heat current. For the 
electronic thermoelectric power eqn. (3) reduces to 
eQ=—c/n sg S ean hal’ (eae oe 
(for a thermodynamic derivation of eqn. (14) see Keesom 1913). We may 
remark here that Joffe (1952) discusses the lattice heat conductivity 
starting from a formula similar to eqn. (12). 

In §3 we give the relevant transport equations. In § 4 we discuss the 
effect of the non-equilibrium of the lattice and in § 5 we give expressions for 
«and @. In § 6 we derive expressions for the electronic contribution to « 
and QY. In §7 and §8 we discuss the influence of the departure from 
equilibrium of the lattice on both the lattice heat current and the electronic 
heat current. Finally, we discuss in § 9 our results and we compare them 
with experiments and with the results obtained by other authors. We 
shall use as far as possible the notation af A or M. 


§3. THe TRANSPORT EQUATION FOR THE ELECTRON DISTRIBUTION 
FUNCTION 


We restrict our discussion to n-type semiconductors and use the model 
introduced by Wilson (1931, see also TM p. 118). It must be remarked, 
however, that our considerations can easily be adapted to the cases of 
hole conduction or of conduction by both holes and electrons. 

We shall assume that the (normalized) electron distribution function 
/(k) differs from its equilibrium value f, and that we may assume it to be 
of the form (M 17, A 1) 


H(k)=fo—(ki/RT) c(n) (Ofo/On). - . . . . « (15) 


We have assumed that there is present an electric field € anda temperature 
gradient d7'/dx both in the x-direction. In eqn. (15) k is the electron 
wave vector and /, its v-component, Rk in Boltzmann’s constant, and 


n= (2—U/RTAY Ta a A 
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where ¢ is the Fermi-level and £ the electron energy for which we use the 
“ quasi-free ’ expression 


H=h*k?/87?m, Ce ee eae awe) 
where h is Planck’s constant, k the absolute magnitude of k and m the 
effective electron mass. For f, we have now instead of eqn. A 5, 

Jo=ez0 (—7). . . . . . . . (18) 


We assume the lattice distribution function N(q) to be given by the 
equation (M 15, A 2) 


N(q)=No—(41/RT)0,(z2)(ON Oz),  . . . . (19) 
where q is the lattice wave ae q, its ee and 
z=hy,/kT, Ce a Pak i tO) 


where hy, is the phonon energy. The index j distinguishes between the 
three possible modes (one longitudinal and two transverse) and the 
equilibrium distribution No is given by the equation 


N =(e'—1)-1 eee hye te ey, 28191) 


The functions c(7) and b,(z) in eqns. (15) and (19) must be determined 
from a set of simultaneous integro-differential equations. The first of 
these equations has the form (compare M 37 or A 8) 


h2 \3/2 EdT O/r ale 
sey (7) ® A Ae | |e er 


ze—* dz 


/O\3 rte 
AM, (7) B*e(n)+ | (Be(m)—eln+2){B+a)} ap > - (22) 
s ~y etl? 
where @ is the Debye temperature, A is given by eqn. M 5a, 
==(477/3)¥® 4 MaRO/3R7C?, . - . . « (23) 
where MV is the mass of an atom in the lattice, a the lattice constant, and C 
a constant giving the interaction between electrons and lattice. The 


quantity M, is inversely proportional to the mean free path J, for impurity 
scattering and is given by the equation (compare TM p. 284) 


M ,=h?/27?mi,. eer ee a sta (24) 

The quantities A and «(z) are given by the equations 
AR EG LOG Od pe ce es) hey eae: Gee 25) 
o.(z)=4tkT2—Dz2*(T/O)?, eer eed (20) 


where D is given by equation A 13, 

D=(62)?h2/16n2ma?. . 2 . . « - + (27) 
Finally, y is given by the equation 

=(2#/D)2(6/T). te 2a - ea ee cee) 
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We may draw attention to a few differences between eqn. (22) and the 
corresponding eqn. A 8 for the case of a metal. First of all, the Fermi 
distribution A 5 is everywhere replaced by the Boltzmann distribution (18). 
Secondly, the limits of integration of the integral on the right hand side 
of eqn. (22) are no longer +0/7', but +y, where y is given by eqn. (28). 
The reasons for this change are discussed by Wilson (TM p. 262, see also 
Frohlich 1936 p. 234). The final difference is between eqns. (24) and 
A 10, but this is only an apparent difference, because eqn. A 10 goes over 
into eqn. (24) by substituting for the residual resistance py the expression 
(see TM p. 285) 

po=3h®/162rme*l.G. % —. ahaa an 


It must be added that the concept of a residual resistance has obviously no 
sense in semiconductor theory since the ordinary resistance increases 
steeply at low temperatures, instead of decreasing as in the case of a 
metal. 

In writing down eqn. (22) we have made the following additional 
assumptions. First of all, we have assumed that we can use the Debye 
approximation for the spectrum of the lattice waves. Secondly, we have 
assumed that we may neglect electron—electron scattering ; electron—- 
boundary scattering can easily be taken into account by suitably changing 
L,. 

In metals /, is a constant which does not depend on the energy. This is, 
however, in most cases not true for semiconductors, since most impurity 
centres will be ionized and thus scattering electrons according to the 
Rutherford scattering formula. If ionized impurity centres are the only 
scattering centres, J, will be given by the equation (Conwell and Weisskopf 
1950) 

l,=(B/rN) [In(1+- Bd?)]-4, -. . . . . . (30) 


where N is the number of impurity centres per unit volume which for an 
impurity semiconductor is temperature dependent and equal to the 
number of carriers, where 2d is the average distance between impurity 
centres, and where B is given by the equation 


B=4x,?H?/e?, 85. ee, es ae Ween es (31) 


with xy the dielectric constant. 

If we want to take into account scattering by non-ionized impurity 
centres, by grain boundaries and by the boundaries of the sample we can 
replace 1, by 1,* given by the equation 


1/l,*= ei Se ae, ey ra 


where /, is the mean free path for scattering by non-ionized impurity 
centres, grain boundaries and sample boundaries which will be energy 
independent. ; ; 

In the following we shall neglect /, and use eqn. (30) for J,. However, 
this may introduce serious errors at low temperatures, since then most 
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impurity centres will be neutral and the number of ionized impurity 
centres will be small. The mean free path J, will then be much larger 
than /) and M, becomes energy independent. 


$4. THe Transport Equation ror THE LatTrice DrsTRrBuTION 
FUNCTION 


Equation (22) gives us one equation for the two unknown functions c() 
and 6,(z). We shall now derive a second equation which is the transport 
equation for the phonon distribution. This equation is (see M 16) 


qidN dT [ON, 
Mog a? dx | O ec’ eed 
where wu, is the lattice wave velocity given by the equation 
Ug=2rv,/q=(47/3)18 a ROA, een er et (34) 


where we have used the Debye approximation which involves the assump- 
tion that uw» is independent of q. In eqn. (33) the index ‘ sc’ refers to the 
change of the lattice distribution function due to scattering ; in many 
papers the index ‘ coll’ is used instead. 

The following scattering mechanisms may play a role (see M p. 485, 
TM p. 292, Frederikse 1953) 


(i) phonon-electron scattering, 
(ii) phonon—boundary scattering (where boundary includes both grain 
boundaries and the boundaries of the sample), 
(iii) phonon—-impurity scattering, and 
(iv) phonon—phonon scattering. 


We have thus for the right hand side of eqn. (33) 
[AN /At] se—[9N /At]er+ [ON /Ot]y+[ON [Ot mp+(ON/Otln. (35) 


We shall briefly discuss the four terms on the right hand side of eqn. (35). 
As far as the phonon-electron scattering is concerned, it is usually assumed 
that only the longitudinal phonons will play a role. Makinson has, 
however, pointed out that this is in general not correct and that it is 
better to assume that transverse and longitudinal waves interact to about 
the same extent with the electrons. We shall adopt this suggestion and 


eae) (en 112 ae oe (30) 
Ni(q)=Na(g)=35 (4) ee 37) 
Ua (2\=b5 (2) =04c) eee ey ee (SS) 


where the C is eqn. (36) is the same as the one in eqn. (23) and where the 
indices ‘L’ and ‘ T’ refer to longitudinal and transverse waves respec- 
tively. It must be added that since it is very difficult to assign a definite 
value to C, the assumption underlying eqn. (36) is unlikely to affect our 
discussion. greatly. | 
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The first term on the right hand side of eqn. (35) satisfies the equation 
(compare Sommerfeld and Bethe 1933, p. 546; M 18) 


ON; 328q,C?2m?a3 iy d E {1 Te Y 
A x 441+ e(n+2) 
| ot lf 9h4Mu,(e?—1) eee 7) an | 2 hq ( 
Joln)—b,6) | 
3273q,C?m2a3 E | ca exp(—) nf [1 ae] 


~ OREM ug(e?—1) L? J) er a: 
x ely-+2)+eln) | exp (—EiK7) 0, (2) J. Se eee 


For the scattering of phonons by boundaries we can introduce a mean 
free path, /), and we have thus for the second term on the right hand side of 


eqn. (35) 


ON. N;—Ny 1% 0,(z) ON, : 
I —— a J = — “ _ ry . . 4( 
| ot i 9 lp kT ly Oz ») 


Since /p is determined by the geometrical configuration of the boundaries, 
we may take it to be independent of q. 

For the scattering of phonons by impurities we can also introduce a 
mean free path, /;, which in this case will not be independent of g, but 
following Peierls (1929) we, have 


11, (gig aidan eae 


where /,° is independent of g, and where do iS the maximum value of ¢ 
given by the equation 
Yo= (3/477) 1/3(277/a) = 27RO/Auy. 2 eel ene cae 


Equation (41) for /; corresponds to the Rayleigh scattering law ; this means 
that it will only be valid as long as 27/q is large compared to the * size ’ of 
the impurities. For large values of q 1; will either be proportional to gq. 
viz., when 27/q is of the order of the impurity size, or be independent of g. 
when 27/q is small compared to the size of ee impurities (see van de 
Hulst 1946). 

Using eqn. (41) we have for the third term on the right hand side of 


eqn. (35) 
a FS eA ia ak W192) (4 40N, 
| at Veoo= Ran! L  ¥OrT (+) So. ht (43) 


We are now left with the last term on the right hand side of eqn. (35). 
We shall use the considerations of Peierls (1929) and Klemens (1951) (see, 
however, Herring 1954 a). From their work it follows that the phonon— 
phonon scattering can also be described by a mean free path, Jp, which is 
independent of g, but which depends on 7’ as follows (Peierls 1929, see 
TM p. 293) 


lp=GO/T, T>O; lp=H exp(0/2T), T<@, . . (44) 
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where G and H may be treated as constants. For a justification of 
eqns. (44) we refer to Klemens’ paper (1951). 

Using eqns. (44) we get for the last term on the right hand side of 
eqn. (35) an equation which is the same as eqn. (40), but with J) replaced 
by lp. As far as the determination of the b,(z) is concerned, we could 
combine Jp and /p in the usual manner, but we shall not do this as their 
temperature dependence is different. 


§5. THE GENERAL EXPRESSIONS FOR THE THERMAL CONDUCTIVITY AND 
THE THERMOELECTRIC POWER 


The electrical current density J, the electronic heat current density w,, 
the lattice heat current density w, and the total heat current density w are 
given by the equations (compare eqns. M 35, M 36, M 24, A 14, and A 157) 


ies eta (He) | eye B3!? ¢(n) de, Ria Se 45) 
Ae a(t) eke B®? ¢(n) ee ee Pha ed (46) 
wy 5 Fy Nia) haul ate | PC) 

ea iat Gah ae 24 b,(2z) ot. (47) 

W=W,+W,g > See ea ot ete 2s (48) 


where we have used eqns. (37), (30), (34) and (42) and the relation 
eA amiel gt bo.” ohn Me cee (49) 


which follows from eqns. (20) and (34). Since in eqn. (47) we are dealing 
with phonons only, the limit of integration is here 0/7 and not y. 
Once again we write (see eqns. M 38, M 37a, A 16) 


2 3 IT Pel 3/2) 
any, (ss:) (7) {Aen + Fe pe semn |}. (0) 


with A given by eqn. (25). Furthermore we write 
Wy = Yd AGO AGL), =. Agee pide. = (D1) 


a result derived in §7. We introduce now the following notation 


Ss ofo 
‘ bids, Bm o(n)22 dn, . . . (52) 
K (m,n) 7 ae. (7) ei 
L;=—7t | | Dee o0(q) Ldn, eet. OS) 
J —C/kT id! 
7=(8A/3rh)(6/T)%, Los tee ee COS) 


+ There is a misprint in eqns. (A 14) and (A 15) in the exponent of 872m/h?. 


P.M. SUPPL,—APRIL 1956 S 


252 D. ter Haar and A. Neaves on the Thermal Conductivity 


and we get from eqns. (45) to (48) and (50) to (54) 


3 3 5 Sh ‘ 
jes K (3,5 )ed+{ & (55 )a ele \ . add eeoe) 
a5 ( ee 
w-(-n(4) am (3D) + [eG8) 2] 
5 if 
4+y E G. 5) pels | +0}5, (56) 
From eqns. (3), (5) and (55) we get for Y the expression 
Q=0.1%, oe Os aoe Toe en 
e(),— —|k @ 5)/ TK (5. ey Lt/T +. eos 
e=—Ly/K (5,5), ee er Ch « (59) 


where the second term on the right-hand side of eqn. (58) derives from the 
fact that A of eqn. (25) differs from the expression in brackets on the 
right-hand side if eqn. (3) by a term —(¢/7')(dT/dx). The splitting up of 
@ into two parts is only done for the sake of convenience. The two 
simultaneous eqns. (22) and (33) are so complicated that one must 
introduce simplifications. What we shall do is the following. We shall 
first of all in §6 calculate the function ¢(7) from eqn. (22) assuming the 
b,(z) to be zero. Then we shall in §7 calculate the functions 6,(z) using for 
the c(7) the expressions derived in §6. Finally we shall use the expressions 
for b,(z) derived in this way to find the c‘”)(7). Afterwards we must 
check whether the expressions for J and w derived in this way are consistent 
with the Kelvin relations, that is, whether for the total expressions (55) 
and (56) satisfy eqn. (10), or whether 


Lg eS a 5 3 ay fe 
K(S. 5) +2L=K (5. 5) —evk (5.5) re (60) 


is satisfied. We shall find that, indeed, this relation is only approximately 
satisfied, but as we have not been able to solve the simultaneous integro- 
differential equations and as the difference between the left-hand side and 
the right-hand side of eqn. (60) is only in the numerical coefficients, we 
felt that our results could probably still be used to indicate the general 
behaviour of the thermal conductivity and the thermoelectric power and 
their dependence on various factors. 

The thermal conductivity « follows from eqns. (45) and (56) in the usual 
manner. The condition J=0 is used to eliminate the electric field ¢ 
and as « is the coefficient of —(d7/dx) in the expression for w we have 


K=KetkytKeg, ty oy ee re ek 

iy the Al 363 W/ 73°. D \ = feos 
«= [K(5. 5) K (5: 5 )-K(5>5)K(5. ai mK (5.3); (62) 
Kyg=—J, SS Re ER a aia) 


-[o 98 {oie 
egg — LK 3°95 [K (Fg . Of ote he” PON wad ter er. Pal otro (64) 
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We have split up the expression for « into three parts which will be 
considered respectively in §§ 6, 7 and 8. 

At this point we want to discuss briefly a mechanism discussed by 
Price (1954, 1955, see also Frohlich and Kittel 1954 and ter Haar 1956) 
which under certain circumstances may give an important contribution 
to the thermal conduction and also to the thermoelectric power. Price 
suggests that there may be an additional term in the heat current due 
to the fact that the carriers may transport ionization energy. The 
contribution «,from this effect to the thermal conductivity is of the form 


Ka=N(R/e)[ ey po/(Mit+ fe) (AYRT), . . . . (65) 


where » is the number density of the carriers, 1, and 1, are the mobilities 
of the holes and electrons, and A is the energy gap, or ionization energy. 
From eqn. (65) it will be seen that «, will be equal to zero as soon as 
one of the mobilities is zero, and hence we should not expect a contri- 
bution from this effect for our model. This can also be seen as follows. 
In terms of the electron distribution function f we have for w and J the 
equations (compare eqns. (45) and (46)) 


J=—eluf de, es ag, Lege oak onan age eae (59 
w=[(4 +4mce?)uf de= A|uf de+4mn| uc?f dc=—(Ad/e)J+w,, (67) 


where u is here the x-component of the electron velocity, ¢ its absolute 
magnitude and dc a volume element of velocity space. Since in order 
to find x, we have the condition J=0, we see from eqn. (67) that in our 
case xa=0. Since neither J nor A are temperature dependent the 
divergence of the first term on the right-hand side of eqn. (67) will be 
zero and from eqns. (7) and (1) it follows that there will be no contribution 
to the thermoelectric power either. 


§ 6. THE ELECTRONIC CONTRIBUTION TO THE THERMAL CONDUCTIVITY 
AND THERMOELECTRIC POWER 


In this section we shall assume 6,(z)=0, or, N(q)=No. This means 
that we are only considering the terms @, and «, in @ and «x. In order 
to evaluate the K(m, n) we must determine the c(y) for which we have 
from eqns. (22) and (50) the following equation 


e+ Ze: * dz 
E"=AM, (jG) Bema) | z | Hon) en +2) (+20) a 
(68) 


We first of all note that the quantity y given by eqn. (28) will practically 
always be small compared to unity. In fact, inserting into eqn. (27) 
numerical values and using m=, a=1A4 we have that y< 1 as long as 


T>3kO?|/D=10-@?~10°x, . . . . . (69) 


$2 
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if we take @—300°K. In deriving inequality (69) we used the fact that 
we are dealing with a Boltzmann distribution for the electrons so that 
we may put H=8rT. ; 

As long as y is smaller than one we can expand the integrand in eqn. (68) 
in powers of z. This is more appropriate since the occurrence of a 
factor e~*| 1—e~*|-!_ means that small values of z will give the main 
contribution to the integral. Physically this corresponds to the fact 
that the main electron—phonon scattering is by low energy, long wave- 
length phonones (compare TM p. 284). 

We now solve eqn. (68) by the method of successive approximations 
(Wilson 1937) and split the problem into two parts according to whether 
the temperature is large or small compared to the Debye temperature. 


(a) High temperatures (O<T) 
In this case we may neglect the first term on the right-hand side of 
eqn. (68) and we write 

o!")(n)=c9'"'(n) Fey '(n)+---, « » « = «9% (70) 
where the expansion is essentially one in powers of 90/7. Inserting 
expression (70) into eqn. (68), expanding the integrand in powers of z 
and retaining only the first non-vanishing term we get, using eqn. (28) 
for y, when we neglect all terms but cy” on the right-hand side of eqn. (70) 
65((n)=(2E"|Dy*)(O/T)?=4DE"-(7/O).. . (71) 
Substituting now the first two terms of the expansion c‘")(z) into 

eqn. (68) we get for c'"’(y) the equation 


T\2 T \2 p+72%e-* dz 
e4(%(n) —=4DE*4(& | | BID 6) Pa eres 


x {1-4 beyr# (1448-25 Jt], $s ee Ls 


where we have used eqn. (71) and where in the integral containing 
€4'"(n) we have made the same’ approximations as the ones leading to 
eqn. (71). In eqn. (72) 6 is given by the equation 


=kT/E, oe Boo SATS ee oak SS ee 
and we see that the following relations hold 
S~w1, <1. nt yaw Laer oe 


In order to evaluate the integral on the right-hand side of eqn. (72) 
we expand the integrand in powers of z. After integrating we are left 
with a power series in y. The first term cancels against the EH inside 
the square bracket on the right-hand side of eqn. (72) and the second 
term gives the result 


€y'")(n) = 3 DE" *92(7T'/O)?| {(—n? + 8n—12)82/12} 
+ {(2n—7)8/24}— (1/18)], (75) 


and we can see that, indeed ¢,'")(7) is small compared to Co\")(n), the ratio 
being of the order of magnitude of y?. 
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eT 


Using eqns. (16) and (18) we have 


~ p 


? Of 
late ci eae) Oxp(G/RE)) 1 (FO) 


and substituting expressions (70), (71) and (75) into eqn. (52) and using 
eqns. (18) and (73) we find for K(m, n) the expression 


‘ 


K(m, n)=3 exp (¢/RT)rD (5) ‘(RT)™+"-2(m-4-n—2) ! 


e) 
2RO? (—n?+8n—12  2n—7 mtn—1)° 
ee Nery 
DT | 12(m+n—2) 24 13 oee.| 


At this point we draw attention to the fact that in contradistinction 
to the case of metals we have now that in general 


K(m, n) AK (n, m). RES Gerd ea el (eS) 


This means that at this point already the Onsager relations, in the form 
of eqn. (60) are violated, and that therefore the numerical coefficients 
of the last term on the right-hand side of eqns. (79) and (84) are not 
reliable. However, one must bear in mind that, anyhow, a theory such 
as the one used in the present paper can never claim to give an accurate 
quantitative picture, as from the beginning one introduces too many 
simplifying assumptions. 
Substituting expression (77) into eqn. (62) we findt 


« =exp ({/RT)7DR°T?(7/O)*[1—(71RO2/72DT)]. . . (79) 


Introducing the carrier density » and the electron—phonon mean free 


path /,,, by the equations (compare ter Haar 1954, p. 246 and TM p. 265) 
N= 2(2rmkT/h?)3’? exp (C/RT), Share tpie (SO) 
Lp=h?DAGO/4r*mT, Sue tes a(S L) 
and using eqn. (54) for 7 we can write eqn. (79) in the form 
Kak ie ROTI TL) ae, emits cee, <a Se) 
where 
Ke == (32/9) 4(nR)(RI mn) th ee te 2. (83) 


is the expression for the thermal conductivity according to the elementary 
electron theory (see Tm p. 232). The right-hand side of eqn. (83) is 
written in such a way that its resemblance to eqn. (12) is immediately 
clear. Since RO<D and since we have been dealing with temperatures 
above the Debye temperature, the correction term in eqn. (82) is always 
very small. 

i es a et 


+The numerical coefficient inside the square bracket is different from the 
one given in a previous communication (ter Haar and Neaves 1954) which 
to) 


was in error. 
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Substituting expression (77) into eqn. (58) we get 
eQ),=(¢/T)—2k+(11k?O7/144DT), . . . . (84) 


where the first two terms on the right-hand side of eqn. (84) are the usual 
expression for the thermoelectric power of a semiconductor (see e.g. 
TM p. 232). The last term is a correction term which also here is always 
small. 

We may draw attention to the fact that for semiconductors neither 
« nor Y are small effects. . 


(b) Low temperatures (T'<@) 

As long as inequality (69) is satisfied we can still evaluate the integral 
on the right-hand side of eqn. (68) by expanding the integrand in power 
of z. For temperatures below the Debye temperature, the first term on 
the right-hand side of eqn. (68) is, however, the dominant one and we 
obtain the zero order approximation to c'”)(y) by neglecting the integral, 


Co (mn) =(E"-2/AM,)(T/O, . . . . . . (83) 


where we have to bear in mind that VM, is not a constant, but depends 
on # (see eqns. (24), (30) and (31)). The first approximation is obtained 
by substituting eqn. (70) in eqn. (68) and neglecting in the integral all 
terms, but c)'")(y). The results is 


cn) = —2(T/0)4(B"-2/A2M 2D). . . 2. . (86) 


In deriving eqn. (86) we have made use of the fact that the term Bd? 
_ under the logarithm sign in eqn. (30) can be written in the form 


Bd? = 4x2E?d2/e#= 4X 1087? ig ?, ss. (87) 


where we have put H=Sk7,, x=10e.s.u.,d=n, V8 (2imp 18 the impurity 
concentration). From eqn. (87) we see that only at the highest impurity 
concentrations (7,,)210'%em~*) and lowest temperatures (7’<5°x) Bd? 
becomes less than one. At these low temperatures, other scattering 
mechanisms will take over, and we may take M, to be constant. 
Otherwise—and this means for practically all cases of interest we may to 
a fair approximation neglect the variation of the logarithmic factor 
with energy, and put under the logarithm H=3k7', This means that 
we may assume I, to be proportional to #~*, at any rate in the approxi- 
mation used here, and we shall write 


M,=M,'(kT/E)?, a eee Re 
L=l,'(H/kT)?, HL a pled Ate ee eee 
where we shall take M,’' and 1,’ to be independent of energy. 
Using eqns. (85), (86) and (88) we get from eqn. (52) 
K(m, n)=(7/AM,')(T/0)3(m-+n) (RT )m+r-2 
exp (€/RT')[1—(m+n—1)(m+n+2)(1,"/lep)], - (90) 
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where we have used eqns. (81), (24), (88 and (89) to write 
Pain carat LO) ee. 6 14 (91) 
Substituting eqn. (90) into eqn. (62) we get + 


ke=(247/AM4’)(T/O)(RT? exp (C/RT)R[1—50(l;'/Iop)] 
= (512/m)¥?(nk)(RT)m)17l,[1—50(L,'/lg]. =... « (92) 


“ 


If 1, had been energy independent, the numerical constant would have 
‘been (32/97)? as in eqn. (83). Since we are at temperatures well below 
the Debye temperature, /,’</,, (compare the neglect of the integral on 
the right-hand side of eqn. (68) and the ratio of c,'")(q) to ¢)”(n) which is 
according to eqns. (85), (86) and (91) /,’/1,,) the correction term in eqn. (92) 
will, in general, be small. As soon as it becomes appreciable we can no 
longer use the approximations applied in this section. 

From eqns. (58) and (90) we get for the thermoelectric power 


Ce) 4h 40K 1,,)0 ot = <3. (98) 


It M, had been energy independent , the term —4k would have been 
—2k as in eqn. (84). While at high temperatures the first two terms 
are of the same order of magnitude, at low temperatures the ¢/7' term will 
dominate. A fortiori, the last’ term on the right-hand side of eqn. (93) 
will always be small. 


$7. THE LaTTice CONTRIBUTION TO THE THERMAL CONDUCTIVITY 
In order to evaluate w, we must know 0,(z) for which we have from 


eqns. (33), (35), (39). (40) and (43) 


4p = 
ba=[p+ p+ png) A+ ee exp cme) |” 


lp 
a 7” ATU 
x | 3ou—« *)| ett rd} (14 ha ) ee Lz) ein) 
_ fu aT (94) 
QnT dx \- 


where we have used eqns. (49) and (42) to express q/q) in terms of z, and 
where @ is given by the equation 


O= 327 C*meRT/ORtMu,?, . . . . - - (95) 
which by using eqns. (34), (23), (81) and (27) can be written in the form 
O= 2/91, Me PPS ee Rees A OO) 


+ The difference between eqn. (92) and the equation given in a previous 
communication (ter Haar and Neaves 1954) is due to the fact that then we 
assumed M, to be energy independent as in the case of metals. 


258 D. ter Haar and A. Neaves on the Thermal Conductivity 
Substituting eqns. (94) and (96) into eqn. (47) we have 


67kO (T\4 O/T ze* dzT 1 ae pf a 
Wy= hat (5) if aan Ete ae iF mG 5) a Tes an é ) 


0 


x exp (cer) | 


1 Liew a 47MUg 
| ane ae inf (14 in Jetn +2)+e(n)} 


hu, dT 
OnE dx 


(97) 


In order to evaluate this expression we must simplify it. We cannot 
neglect the term with the c(j), but we shall assume that in different 
ranges of temperature one of the scattering mechanisms will dominate 
and that we may neglect all but the dominant scattering mechanism in 
each of these temperature ranges. In order of increasing temperature 
the dominant scattering mechanisms will be (see M) : phonon—boundary 
scattering (range I), phonon-electron scattering (range I1), phonon— 
impurity scattering (range III), and phonon—phonon scattering (range IV). 
Although the Debye temperature will probably occur somewhere between 
range II and range IV (see TM p. 294) we shall for ranges II and IIT 
consider two cases: Ila and IIla, assuming the Debye temperature to 
be higher than the temperatures in those ranges, and IIb and IIIb, 
assuming the Debye temperature to be lower than the temperatures in 
these ranges. We shall in the following, with Makinson, neglect the 
term with 47mu,)/hq and we refer to M for a justification. 


(i) Temperature range 1 


Since we are at temperatures below the Debye temperature we can use 
for the c(7) the solution given by eqn. (85). We neglect here ¢,!"(n) 
since this only amounts to a small correction as we saw in § 6. 

We first of all evaluate the integral over 7 in eqn. (97) and we replace 
¢(y+2) by ¢(7) which is allowed since the main contribution to the integral 
over z will come from small values of z. The result is 


rad a oe 
e—"dy | (: fe wt) o(+2)-+e(n) | — 


es 
Ss h2 3/2 -00 1 dT (2 
se 33/2 pn 5/2 9 
Myke (2) 4 | — RT 5 ant ip T dx oe ; "an} 


h2 1/2 1 
= —2), (<aner) exp (RD) 4 (5) a+er(; rag . (98) 


where r(3)= St? and r(; 


bol ws 


) =1571/?/8 (e.g. Larsen 1948, p. 251), 
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Substituting expression (98) into eqn. (97) and neglecting in the 
denominator all terms but 1/l, we get after integration over z for w, 
the equation 


“ etn (Lat he (RT Vat ales dT 
7 ha \O | ( a) 12lep per. (RT) {249k | 16/7) 
hu k dT 
eee, aes! oe tn 5 eee et 05 (99) 
where we have used eqns. (24), (88) and (89) and where 
Olp 
LAO) | ee Mai dz0 2 eS... -(100) 
J0 
\ rOlT 
J(O/L)= | ere(e = 1ystde. ws... (101) 
J0 


These factors will be approximated by us as follows (see TM p. 336) 
I, (O/T )=1,,(c0)=n! C(n+1), J,(O/T)=J,,(00)=n! C(n), O>T, . (102) 
LAO =a (O 1)" I (Oj =(n— lA (O/T), O<T, --,.- (103) 


where C(x) is Riemann’s zeta-function (C(v)=2',n“*). 

We can use eqn. (55) to eliminate A in terms of J and d7'/dx and using 
eqn. (80) for n and (42) for @ we find for the coefficients y,; and g ;(=—x«,, 
see eqn. (63)) 


yy= —nbyl, RTT ,/121,,m2u? eK (2, 3)... (104) 
Kg=—G1= (31, RUT ,/a)(T/@,  . . . . (105) 


where in the last equation we have neglected a term of the order /,’//., 
as compared to the term retained. 

We can rewrite the last equation by introducing the lattice specific 
heat c, by the equation (compare Sommerfeld and Bethe 1933, p. 547, 
we have an extra factor 3 to be taken into account the fact that we assume 
the two transverse modes and the longitudinal mode to be equivalent) 


cy = (3/20) da )E_N (q)hugg=(3huy/160°7) E 2eq3(e—1)- dq 
ee (OR CLO) kasi te. 9-1). ea ee et ne ome em LOG) 
where we have used eqns. (42). Combining eqns. (105) and (106) we have 
EZ ECU liege Nee on ot elMotiter + saul Ot) 
and similarly using eqns. (42), (54), (80), (90) and (106) we find 
Yy=—(7eg1, TL ,/144nl,,eJ,) exp (C/RT). . . ~ (108) 
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(ii) Temperature range Ila 


We shall still use the result of eqn. (98) and now neglecting all terms 
in the denominator but the one involving J, 1 we get for wy. 


27mlep h2kRT 1/2 Ll. , ahh @Q 
 Cpeeethaty (5) | (ar 12ep pce V4 2D 


+e exp (tik) Ka(z) |) Meee 
where . rOlT VS Pe 
K,(6/T)= |. 2me-*(1—e*) 8 dz. =. . . . (110) 
Eliminating once again A we find now 
Yrs=— TCT /32ne;, we ae ee 
Kg= —Jtta= Cg Uolep CXP (—C/RT) (Ka/J4), . . ~ (112) 


when once again we have neglected a term of the order /,’/1,, in « and 
where in ¥;;, we have again simplified the final expression. 


(iii) Temperature range Illa 


We now retain only the term with (/;)~1 in the denominator, and 
we get for w,, Yi, and K,, 


67h O ( (h?RT\12 ote as () 


hu, kdT (0) 
= ae dx — J CM . . ° . * ° ° . . . . . (113) 
Yara = — (meg TT 9/1440, eJ ,)(9/T)* exp (C/RT), . . (114) 
= —Jrtta= FC gU gly.” (O/T )*(J o/J 4). Ce Rey Wet ncn (115) 


(iv) L'emperature range IIb 


We now must use for the c(” (7) the cy‘")(7) given by eqn. (71). Using the 
same approximations as the ones leading to eqn. (98) we have now 


ie ¥ 4mu h2 \12 
ui d ‘0 > hat —— 
— C/kT 1 | (1 Tr ae ee) che aa 2lep (<aen) 


exp (¢/RT7') ri) A+kI(3 a ee tant LLG) 


where I'($)=7/2, ['(3)—71/2/2 (e.g. Larsen 1948, p. 251). 


We can use eqn. (116) and neglecting in the dominator of. (97) all terms, 
but the one ee ly } we obtain wy, Yr and Ky, 


ae hekT\12 iT\ , (9 
Wy= ha? 6) la(sa) {24-4 daz ) 


hu,k At (0) 

+> exp (—{/kT) —— -K.(z) |, Sve ek Gem ITT} 
Yi»= —7¢, T'/4ne, b Welw, 0 AKU Pe eee ea Dae 
K9=—IJt1p= 3, Uolep EXP (—L/RT)(Ky/(Ja), - = - (119) 
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where now we have used eqn. (77) for the K(m,n) and where in kK, we have 
neglected a term of the order of magnitude (h?/67mkO)?(n/a)(0/T) which 
for n<10'8 cm-*, T>@=300°K, a=3A4, m=m, is smaller than 0-08. 


{v) Temperature range IIIb 


We now retain only the term with (J,)-1 in the denominator and with 
the same approximations as in the previous subsections we get for wz, 


Yup and key 


6715 OT 1 (h2RT\122 ' dT C) 


8am 
hu, k dT 6) 
4+ om) | ee Mere: 0-8 -T77(1 20) 
Yrrp= — (76,17 T/18n el, J 4) (O/T)4 exp (C/RT), .. . (121) 
Kg —GJrrp= 36g Mol (O/T) (Jo/J 4). eee Pee (22) 


(vi) Temperature range IV 


Finally, we retain the term with /p~1 in the denominator and we now 
get for w,, yyy and «,, 


67lp (T\3 71 (h?RT\12 IT’ @ 
w= — Fe? (B) Ee exp (RT) {24 +R F }La() 


_ huykdT C) 
Yry= — (me glp TI ,/18nel.,J 4) exp (C/RT), . . . (124) 
Ky= —Gry=FCgUolp- 5 : . . . . ° ‘ 5 (125) 


§ 8. THE CHANGE IN THE CURRENT DENSITIES DUE TO LATTICE 
NON-EQUILIBRIUM 


We must now evaluate the contributions to « and QY due to ¢”)(7). 
We shall assume, as we mentioned earlier, that in evaluating these 
functions we may use for the b,(z) the expressions obtained by assuming 
c\)(n) to be zero. We appreciate fully that we should really obtain 
solutions of the simultaneous integro-differential equations, but we have 
not succeeded either in doing this or in estimating the errors made in 
our calculations. However, for reasons to be gone into at the end of 
this section we feel that our results are probably not too bad and can be 
used for comparison with experimental data. We shall again distinguish 
between the temperature ranges I, Ha, IIa, IIb, IIb, 1V. In the first 
three cases, the ratio of the first term in curly brackets on the right-hand 
side of eqn. (94) to the second term is of the order of J,’//,, and thus 
negligibly small, while in the last three cases its ratio is less than 0-1 
for n<1018 em-3 (compare the remarks at the end of subsection (iv) of § 7). 
We shall thus use the following expression for b,(z) in the present section, 

hugal li. 71 OS: iets se pee iB 
sel hk ie an(g) @ Te e *) exp (C/RT) : 

(126) 


b,(z)= 
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We now introduce a quantity B(n) by the equation (compare eqns. 
A27 and M37a) 


Oe 2 \3/2 3 
Preyat? ae i po 4(san) (7) Bin) cyte. 
so that c'”)(y) is the solution of the equation 
2 bi 2%e-* dz 
AU eek () BPc\(m) + [EH )—cl( +2) {EB +a(2) @) | a 
(128) 


For temperatures above the Debye temperatures we may neglect the 
first term on the right-hand side of eqn. (128) and we get after the usual 
expansion in powers of z the integral sign and the neglecting of all but the 
dominant term (compare the derivation of eqn. ‘ 1)) 


cn) =1D(T/O)2B-2B(n ee bt 


For temperatures below the Debye ei the integral can be 
neglected and we get, using eqn. (91), 


(yn) = (AM, B?)(1/0)°B(n) = D(T/O)( RL) *(l,'/lep) )B(n)-  - — (130) 


We can now consider the six temperature ranges. We use eqn. (64) 
for «,,, (59) for Q, and (53) for the L;. In the temperature ranges I, La, 
IIIa we use eqn. (90) for the K(m,n) and (130) for c(n), and in the 
other three temperature ranges we use eqn. (77) for the AK(m, n) and 
(129) for c‘”(y), in both cases neglecting the correction terms in the 
expressions for the A(m,n). Finally we only retain the term with 
ly + in the denominator in eqn. (126) in temperature range I, only the 
term with Jey in ranges Ila and IIb, only the term with (J;)-! in 
ranges IIIa and IIIb, and only the term with /)~! in range IV. We 
then get for B(m), c\ (7), L;, «,, and Q, in these six cases the following 
expressions, where in eqn. (127) we have neglected the first term in «(z) 
in accordance with inequality (69). 


‘ 
(i) Temperature range I 


hu, ly D (872m\3/2 
B= ear (ae) > (131) 
where we have neglected terms of the order (T'/@)°. 
hu ly D? ha S72m\3/2 
ol) 0 °bD By: 
= Sn TAR © Toy (Fe ) - (132) 


c,ly ( 2 \12/0\5 Sy Se ae eee 


ep 


cy ly (2akT\12 1,’ (0\5 
ee | ) (5) exp (C/RT) 0s aoe) us a tae) 


mM 


20s TC, ty h2 \3/2 /@\5 
T1920, lep \Q2armkT my FPSO a, Be ie see Lae 
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(u) Temperature range Ila 


9R Uy lep D (87?m\3/2 (T\5 
B(n)= Terr (ae) (5) exp (-—C/RT) (2J,+],), 


0) 9h vg Ll, D2 (82?m\3/2/7 
OND = Bin @ ARO \ RE ( 


a 


ihe So TE , 
be 25 (ep) Fab L OP C+B21RT I 


9, ly! (2akT\12 
et 2 eG Ay \ onan, ; 
Deo) ome 
ea | 2-2 - 
0, | +3 | 128n 


(a) Temperature range Illa 
h Uo i; D [ 872m \3/2 
P= STA Ae ) ; 
— EE a 872m \3/2 
= Ten TAR Oleg \ re) ° 


7 €, °°) 2 \12/@\51,’ SEC a ae 
ae s( ser) (7) ie exp (C/RT) i +-3/2)(kT) 


cl (IckT\121 (O\5 
af oo (=) -(7) COD 


! re h2 3/21; (/@\5 
©0,=— say. sener) ia (P) 


(iv) Temperature range IIb 


9A Uy ley D ( T \5 ( 80?m \3/2__ 
B= (B) (Fr exp (—C/RT') J, 


where we have neglected a term of the order of 0/7’, 


2 a 2 3/2 
ing telen De (2 (S28 7 oxy (tn, 


167 T AEH?\ © h? 


L j= 5 Cglep(}mkT)°T(-—J)(RT) -, 


3 
Keg —_— 3 ¢,(7RT/2m)lep, 


37 Cy 
Wg 16n ° 


5) es? (—C/kRT) (24 ,+-1;,), 
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(146) 
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(v) Temperature range IIIb 


hu, © D ( 8x?m 
a= i (Far ) to 
hu, 1; D? Srr2m)\ 3/2 
on) = roms) (5 fe ) it 
(0) V2/@\41 a 
t= (ea) (gp) perp CRE) HUNT) 
at! 
c, 4 (akT mn giex 


ic T Cy 1, 3/2 7) 7 
ed = — 48 lep — 


(vi) Temperature range IV 


hug lp D [ 82?m \3/2 (T \5 I 
= aTA\ re) \o) 7# 
hu, lp D2? (T \7 ( 822m \3/2 
n= ertrm g) (oar). te 


© keg tf Ses ee 
sigs (<en) sex (a2) PERT), 


Bee) 


C, lp (m7RT\ZT, ; 
Sere 2 12 ( oa) 7, oP (C/RT) 


vin Cc; lp h? 3/2 I, 
es aR Te (ser) ah 


(159) 


(160) 


Before we start the discussion of our results, we want to report briefly 
the results obtained for the case of metals, and also we wish to see whether 
or not eqn. (60) is approximately satisfied. This equation is not 
satisfied in the case of metals, since to the approximation used in our 


present paper y=0 for metals, while L;=0 and K(m, n)=K(n, m). 


We 


give therefore our results for metals including only an approximate 


value for the numerical coefficients involved. Our results were 


eQ,= —0:3(c,/n)(Lp/lep), cle Gas ok th: os A rene 
eQ,=—0-3c,/n, : (range Ila) 
eQ) = —0-002(c,/n)(1; (Lep)(O/T')4, (range IIIa) 
C6) EDC IM, Ao ass) a ek aes a ene (range IIb) 
e),=—0-1(c,/n)(ly (lep)(O/7), (range IIIb) 
e(), = —0-05(c,/n)(Lp/lep)(0/7). (range IV) 


(161) 
(162) 
(163) 
(164) 
(165) 
(166) 
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We may draw attention to the fact that eqns. (162) and (164) are, 
indeed, similar to the expressions derived in A, apart from the numerical 
factor and the sign. The different numerical factor arises from using 
Wilson’s method of integration instead of the approximate one used in A 
and the different sign arises from our use of the Thomson coefficient 
instead of @ (see Sondheimer 1956). The results are not sufficiently 
reliable for comparison with experimental data. As far as eqn. (60) is 
concerned, to a fair approximation K(m,n)—K(n, m) as we saw in §6 
so that we must check whether or not - 

ple CU K (Ge) 06 QO) y/T 8.) a.e. (167) 
Comparing the results of the present section with those of § 7, we see that 
eqn. (167) is, indeed, to a first approximation satisfied, apart from the 
occurrence of factors (0/7) in the low temperature cases I and IIIa. 
There are indications that these factors should not be present in a more 
satisfactory theory and we have therefore discarded them for the discussion 
in §9- 

§ 9. Discussion 

We can now put all our results together using eqn. (61) for « and 
eqn. (57) for Y. We first of all notice that the expressions for « derived 
in §8 are all of the same order of magnitude as the terms neglected in 
§7, that is, they are small compared to the main contribution to x. 
We shall therefore neglect «,, in our discussion. As far as x, and Q, 
are concerned, the correction terms are small compared to the main 
term in all cases of interest so that we shall also neglect those in the 
present section. We can now combine the results of §§ 6 to 9 as follows. 


K=10(nk)(RT/m)4/1,’++0-3¢,uglp, - - . - - (range I) (168) 
_K=10(nk)(RT/m)*71,'+-1-5¢,uol., exp (—C/RT),  . (range IIa) (169) 
ied (RT /m)21, -0-Olesugl OP, - +. (range IIIa) (170) 
= (nk)(RT/m)*/?1,, + 2-3¢ Ugly exp (—E/RT) 0), (range IIb) (171) 
K=(nk)(RT/m)*lep+0-3¢,u hi, . . = . (range IIIb) (172) 
K=(nk)(RT/m)*lep+0-3¢,Uglp, -. - : (range IV) (173) 
eQ=(6/T)—4k—0-001(c,/n)(Ip/lep) exp (ik), (range I) (174) 
eQ=(C/T)—4k—0-2(c,/n), . . . (range IIa) (175) 
ae (¢/T)—4k—0-001(c,/n)(L (i: |lep) € exp GD), (range IIIa) (176) 
=(¢/T)— 2k—0-6(c,/n) Sane . . (range IIb) . (177) 
—(¢/T')—2k—0-6(c,/n) (1,/lep) exp (t/kT),. Rh (range IIIb) (178) 
=(£/T)—2k—0-1(c,/n)(Lp/lep) exp (S/RT)(O/T), . (range IV) (179) 


oe we have only given approximate values of the numerical coefficients 
and where in eqns. (174) and (176) we have dropped the extra factors 
@/T in accordance with what we said at the end of §8. In eqns. (170) 
and (176) /,@” is the average phonon-impurity mean free path given by 


the equation 
La (O/T jt. sets ok Seen eee (ESO) 
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First of all, we shall study the general temperature dependence of th , 
various expressions. Assuming the number of impurities to be propor- 
tional to n, 1,’ is proportional to 7”, while Jep is proportional to T-! {606 
eqns. (89) and (81)). As long as we are not at too low temperatures, 
n and ¢/T' are practically constant. It follows that below the Debye 
temperature x, is proportional to 7°’? and above the Debye temperature 
x, is proportional to 7-1”, while Q, varies only slowly with temperature 
(for its temperature variation we can refer to the papers by Johnson and 
Lark-Horovitz (1953), Lautz (1953) and Fukuroi and Tanuma (1952) 
and also a forthcoming survey article by Johnson). To study the 
temperature dependence of the lattice contributions we take Jp, 1 and 
u, to be constant, Jp) and Jep to be proportional to 7-1, c to be propor- 
tional to 7? below @ and to be constant above 9, n and ¢/T' to be connected 
through eqn. (80), and » to be approximately constant. If we write 
k, proportional to 7” and @, proportional to 7° we get for r and s in the 
six temperature ranges under consideration 


. . once * aaa, ie . Ss . See . Q 
y=3,  ye=7/2, Trmo=—1, Tp =3/2, Tum=9,  tv=—1l; . (182) 
ab sath i cy ey ee ee eee > 
$;=5/2, 8y_=3,  Sqya=— 3/2, Sqp=9, Ss p=— 2, Sty=—4/2. . (182) 


From eqns. (181) and (182) we see that both x, and Q, should show a 
maximum and this has, indeed, been observed (see Rosenberg 1954, 
1955, Estermann and Zimmerman 1952, Goff and Keesom 1954, and 
Frederikse 1953 for «, [compare also Goldsmid 1956] and Frederikse 1953, 
Geballe and Hull 1954, 1955, Mooser and Woods 1955, Fukuroi ef al. 
1949 for Q). It is of interest to note that Rosenberg in a very high 
purity sample of germanium finds, indeed, an exponent larger than the 
usual 5/2, although the value of 4-4 found by him is much larger than any 
of the exponents mentioned in eqn. (181). Another point of interest in 
this connection is the dependence of « on n: it is found to decrease with 
increasing number of carriers in accordance with eqn. (169). The region 
where the dependence on » occurs is the same as the region where the 
maximun of @, seems to occur and may thus well correspond to our 
temperature range Ila. This point needs further investigation, however, 
as in the case studied x is so small as to practically exclude the possibility 
of phonon-electron scattering. We may refer to a paper by Ziman 
(1956). As far as the order of magnitude of «, is concerned, using 
O=362°K (Keesom and Pearlman 1953) and a=3A, we find that at 20% 
for n=10'% cm-3lep) should be 10-> cm which does not seem to be un- 
reasonable. At 20°K we find that x,/«, is of the order of 10 forn=10!6 em-3 
and of the order of 1000 for n=1015 em-3. 

As to the thermoelectric power, we also find there an increase with 
decreasing ” as is to be expected. The qualitative agreement with the 
experimental data is quite satisfactory, but quantitatively, eqn. (175) 
gives values for @, which are too large by about two orders of magnitude, 
At temperatures below 0, Q, dominates, but above 8, Q, is the larger, 
which is satisfactory since, as we mentioned earlier, at high temperatures 
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the theory of @, can satisfactorily explain the observed behaviour of the 
thermoelectric power. 

In conclusion we should like to discuss briefly the relation between the 
theory developed in the present paper and those of other authors. The 
theories of Klemens (1954) Parrott (1954) and MacDonald (1954) are 
qualitative and more closely related to a kinetic discussion (ter Haar 
1956). Gurevich (1945) only deals with metals and with @,. He con- 
siders four temperature regions, two below and two above the Debye 
temperature and finds the exponent s to have the values 1, 3, 1, —1 
respectively while we find 4, 3, 0, —1 in rough agreement with his results, 
but not exactly the same. Our results are in direct contradiction to the 
opinion expressed by Pikus (1951) who claims that Q, will always be much 
smaller than Q,. As his conclusion is based on some unpublished 
calculations according to which J.p>T/p it is difficult to see where Pikus’ 
calculations go wrong, since they must have gone wrong somewhere as the 
experimental data indicate Jep</lp) at low temperatures, while Pikus 
claims /,,,>lp for T>3°x. 

Finally, we must say a few words about Herring’s theory (1954). 
In his very elaborate analysis Herring uses what he calls the [/-approach, 
where one calculates the Peltier heat rather than the Seebeck effect. 
It seems to us that his approach does not satisfy the Onsager relations, 
as he assumes implicitly (e.g. on p. 1170) that one may use the equilibrium 
lattice distribution function in calculating the non-equilibrium carrier 
distribution function. He uses the quasi-free model, but discusses in 
an appendix some of the consequences of the many-valley model. It is 
difficult to make a closer comparison between his and our theories as he, 
of necessity, introduces a number of simplifying assumptions which are 
partly different from the simplifications introduced by us. One of his 
results, however, seems to be to us incorrect, namely, that the increase 
of Q, at about 20°K would not depend on n. This seems to us to be both 
in disagreement with Frederikse’s experimental results and with the result 
obtained from a kinetic argument and we feel that our result (175) in 
the neighbourhood of the maximum—at any rate in so far as its depend- 
ence on 7 is concerned—should be more reliable. His analysis goes much 
further than ours in also considering degenerate semiconductors, but our 
analysis is more elaborate, we feel, in its discussion of the various 
scattering mechanisms. 

After this paper was written, Dr. V. A. Johnson of Purdue University 
kindly sent us a preprint of a review article on the Seebeck effect in 
semiconductors which will be published shortly. Her approach is more 
closely related to that by Herring than to ours ; moreover, she is mainly 
concerned with Q,. 
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D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. I. 


Oriented hexagonal-based crystallites of sodium chloride upon mica. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, 


Ple2: 


Oriented square-based crystallites of potassium bromide upon mica, 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 3. 


Oriented silver bromide upon a silver (111) surface. Asymmetrical intensitiy 
distribution about central line indicates pronounced single positioning. 
Ag [110] azimuth. 


Fig. 13 


Pattern from a (100) film of silver chloride showing the formation of 
oriented silver during electron bombardment. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, PI. 4. 


A layer of about 24 of lead on an atomically flat (111) silver surface. 
Ag [110] azimuth. 


Fig. 18 


A layer of about 2A of copper on an atomically flat (111) silver surface 
Ag [110] azimuth, Diffuse copper spots appear just outside vertical 
streaks due to silver. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 5. 


Fig. 19 


A layer of about 24 of silver bromide on an atomically flat (111) silver surface. 
Ag [110] azimuth. Spots due to bromide ; vertical streaks due to silver 
and to secondary diffraction. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 6, 


Fig. 21 

acre 
DAS 
b* 

7 

v. 
Ao? 

Vv 
v é 


s 


pe? oh ae 46 uPA <a> : 


x 
qe 


Oriented triangular-based pyramids of ammonium iodide upon mica, 


DNV A PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 7, 
Fig. 22 


i le Kikuchi line pattern from an oriented silver layer on mica. 
Sagi Se [100] azimuth. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 8. 


Silver bromide layer on a (110) silver surface. AgBr [110] azimuth. 


Fig. 24 


As for fig. 23 except AgBr [001] azimuth, 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 9. 


A layer of about 2 of silver on a heated rocksalt surface. NaCl [011] azimuth. 


D. W. PASHLEY Phil. Mag. Suppl. Vol. 5, No. 18, Pl. 10. 
Fig. 30 


As for fig. 29 except silver thickness of about 20 A, 


Fig. 31 


As for fig. 30 except silver thickness of about 1504. 
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